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Abstract 

a 

We consider the possibility that along the thermal history of the Universe, dark matter 
' (DM) would have been created from Standard Model particles, either through a kinetic mixing 

portal to an extra U{1) gauge field, or through the Higgs portal. Depending solely on the DM 
particle mass, on the portal and on the DM hidden sector interaction, we show how the observed 
I— I DM relic density can be obtained. There are four possible freeze-in/reannihilation/freeze-out 

regimes, which together result in a simple characteristic relic density phase diagram, with the 
shape of a "Mesa". In the case of the kinetic mixing portal, we show that, unlike other freeze-in 
^ I scenarios discussed in the literature, the freeze-in regime can be probed by forthcoming DM 

O direct detection experiments. These results are well representative of any scenario where a DM 

I hidden sector would be created out of the Standard Model sector. 

> 
On 
p 

CN 1 Introduction 

^~~] There is an all series of gravitational evidences for the existence of dark matter (DM). In particular, 

^ from observations of the Cosmic Microwave Background anisotropies and large-scale structures, it 

• ^ appears that a large fraction of the energy density of the Universe is made of DM, with ^dm = 

^ 0.229 lb 0.015 [T]. To explain such a proportion, the WIMP candidates are particularly interesting 

because, by producing thermally the DM from the thermal bath of the Universe, through DM 
pair creation and annihilation, they lead to a relic density which depends only on the masses and 
interactions involved at the time of the decoupling of these processes. Moreover it turns out that for 
interactions of order the electroweak ones, or more generally for couplings of order unity, they lead 
to the observed relic density for a DM mass of the order of the electroweak scale. This offers the 
possibility of DM production at colliders, on top of the possibilities of direct and indirect detection. 

Another simple mechanism to produce the DM relic density, which was pushed forward recently 
O [31 m |5] , is the freeze- in mechanism. Here the relic density, instead of resulting from the thermal 
freeze-out of the DM annihilation process, stems solely from the DM creation processes, for instance 
through annihilation into DM particles, AA — )• DMDM, or through a decay process A — )• DMB. 
These processes, that are assumed to be out of thermal equilibrium, freeze when the temperature 
T drops below the mass of the source particle A or DM particle, i.e. when the production rate gets 
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Boltzmann suppressed. While being out-of-equilibrium, this mechanism is actually also thermal in 
the sense that the source particle A is assumed to be in thermal equilibrium, so that, here too, the 
number of particles produced depends only on the masses and couplings involved in the DM creation 
process. A difference though with the usual freeze-out scenario is that the relic density may also 
receive a contribution from a "primordial" density, which is not washed-out by any process since the 
DM creation process is out of thermal equilibrium. However there exists a whole regime for which 
the DM density produced from the Standard Model (SM) sector can be the dominant one 0This IS 
the possibility we will consider in the present work. Another important difference with freeze-out 
is that the characteristic coupling required to produce the relic density through freeze-in is tiny, 
typically about 8-10 orders of magnitude below that required by freeze-out, implying that the DM 
basically lies in a hidden sector that is feebly coupled to the SM one. Hence it is a priori very 
difficult to probe experimentally the freeze-in scenario, either at colliders or from direct /indirect 
detection. There is nevertheless one way one could consider to test this mechanism in specific 
scenarios [Sj HI [5]: if freeze-in proceeds through A — t- DMB decay (assuming typically that A and 
DM particles are both odd under a Z2 symmetry so that DM is stable) and if, in the visible sector, 
A has substantial coupling to the SM particles, then A particles may be abundantly produced at 
colliders. Their slow decay to DM from the feeble interaction could then be probed from the search 
of long-lived particles at colliders. To be efficiently probed, this scenario requires that B itself is 
substantially coupled to SM particles (being typically charged or colored), so that only DM particles 
are feebly coupled to all other particles. 

In the following we show that there is another possibility of testing a DM relic density created 
through freeze-in, i.e. direct detection. The scenario we consider is quite different. In particular 
we do not assume any new particles in the visible sector, but the ones from the SM. This means 
that in our framework the A particle is a SM particle, which implies that the predicted DM relic 
density depends neither on the mass of any new hypothetical particles in the visible sector, nor on 
the interactions these particles may have with the rest of the visible sector. The DM relic density 
depends only on the DM particle mass, on the connector interaction and on any interaction which 
may be relevant within the hidden sector. One class of models we consider is based on the possible 
existence of a light mediator particle. We will show that such a light mediator interaction, as tiny as 
the one required for DM production through freeze-in, could be probed by direct detection, because 
the elastic cross section with nuclei, which exhibits a collinear divergence, is strongly enhanced. 

The light mediator DM model we consider in the following as an example is very simple. It 
involves, in addition to the SM (the visible sector), a hidden sector composed of a single particle, 
which is assumed to be charged under a new C/(l)' gauge symmetry. This hidden sector is coupled to 
the SM sector through kinetic mixing, C 3 -leFi^^F'^^ (see Ref. [6l|7] for a "secluded" DM structure 
of this type with a massive 7', see also the related model in Ref. [8]). Thus DM may interact 
with the SM particles through a massless mediator, or a light mediator if the C/(l)' gauge group 
is slightly broken (either through the Higgs or Stiickelberg mechanisms). On top of the connector 
interaction, this model involves QED-like interactions between the 7' and the charged DM particle in 
the hidden sector. Another possible motivation for this model is that there are not that many simple 
mechanisms available to stabilize DM (see e.g. [9]). Among these, one of the simplest possibilities is 
to consider the lightest particle charged under an extra, unbroken gauge symmetry. Hence the model 
we invoke here may also be seen as the simplest and yet generic representative of a larger class of 
models. In particular we expect that the features we will discuss should hold for many models of the 
same kind. A further motivation is that, although the cosmological and astrophysical consequences 
of such a new long range gauge interaction have been considered extensively in Refs. |101 [TT| 112] . 

^If both sectors are feebly connected, one could argue that the inflaton is likely to belong to only one of the sector, 
so that it couple dominantly to one sector and reheating occurs dominantly in this sector. 
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this has been done basically without considering the possibility of kinetic mixing connecting both 
sectors. 

Beyond the fact that the freeze-in mechanism may be testable within this specific scenario, the 
main, more general purpose of our work is to determine the various ways a viable DM relic density 
can be obtained in the context of particle creation through a portal. To this end we have explored 
the full parameter space of this model. If both the kinetic mixing and the hidden sector interactions 
are feeble, the DM relic density can be calculated by simply counting the number of DM particle 
produced through the connector. This is precisely the freeze-in scenario. However, if instead either 
one or both of the interactions are stronger, the DM relic density is affected by the thermalization 
of the corresponding processes. We show that this leads to a simple characteristic "Mesa" phase 
diagram for the relic density as a function of the connector and hidden sector couplings. As we 
explain in details this phase diagram displays four distinct phases: (i) freeze-in, (ii) simultaneous 
freezing of the (in-equilibrium) hidden sector interactions and the (out-of-equilibrium) connector, 
i.e. reannihilation, (iii) freeze-out of the hidden sector interaction, and (iv) freeze-out of connector 
interaction. An interesting feature of this diagram is that, in each phase, the relic density is 
essentially characterized by only one interaction, either the connector or that in the hidden sector. 

In order to address the generality of this classification, we consider the creation of DM through 
another simple portal between the SM and a hidden DM sectors, the so-called Higgs portal. We 
show that DM production from SM particles through the Higgs portal leads to a similar simple 
characteristic phase diagram which also has characteristic "Mesa" shape. This holds even though 
the mediator, i.e. the Higgs boson, is massive and may even be heavier than the DM particle, unlike 
for the kinetic mixing where we consider the opposite limit of a mediator (much) lighter than the 
DM candidate. In this sense the "Mesa" diagram is a generic feature of models within which DM is 
created out of SM particles through a portal (or more generally from particles that are in thermal 
equilibrium with the SM sector when they create the DM particles). 

The plan of this article is as follows. In Section 2 we first compute the transfer of energy from 
the visible sector to the hidden sector in the case of the kinetic mixing portal. This is necessary to 
compute the DM relic density in all regimes, except in the freeze-in one. We then compute the DM 
number density produced. In Section 3 we present and discuss in details the corresponding phase 
diagram obtained and give analytical approximations for the relic density. In Section 4 we discuss 
the possibility of direct detection signal through kinetic mixing, and the prospect for testing the 
freeze-in mechanism. In Section 5 we discuss the compatibility of the relic abundance and direct 
detection signatures obtained with the various existing cosmological and astrophysical constraints 
on DM interacting through long range interactions. In Section 6 we present our results for the 
Higgs portal, and discuss in details the similarities and differences between the cases of massive and 
massless mediators. Finally, we summarize our results in the last section. 

2 DM number density and energy transfer from the visible to the 
hidden sector through kinetic mixing 

It is remarkable that the addition of the simplest gauge structure one can think of, that is to say of 
a U{iy gauge structure with a single charged particle, leads to a viable DM candidate |10 | [T H [T2|. 
In the following we will consider a charged fermion (the DM candidate), 

C = CsM + ip'ii^' -m^)ip' , (1) 

with D'^ = + ie'A'^, but the results would be essentially the same for a scalar particle (except 
in this case for a possible additional Higgs portal interaction between 2 DM particles and 2 Higgs 
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doublets, which will be considered in Section 6). This QED' model contains, on top of the SM 
particles, an extra massless gauge boson 7' (the hidden photon), and an extra fermion, that we call 
"e"' (singlet of the SM gauge group), and nothing else. This set up can be coupled to the SM in an 
unique way through the kinetic mixing portal |13| [H] 

^3-|i^r^;.- (2) 

where Fy'^ is the hypercharge field strength. Consequently this model Lagrangian involves only 
three new parameters: a' = e'^/47r, e and m-e' = rriDM- The non-canonical kinetic term implies a 
mixing between the visible and hidden sector photons, 7 and 7', as well as between the Z boson 
and 7'. For all purposes it is convenient to work in a basis where the kinetic terms are canonical 
which can be obtained from a non-unitary transformation. We use 

(3) 
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where 9w is the Weinberg angle and the 9^ mixing angle satisfies tan 9^ = tan 9w / \/l — e^. 
We then have the following couplings to the SM and DM currents 
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and with Jgj^j and J^j^ the corresponding U{l)em and 11(1)' currents. 

As the visible and hidden photons are degenerate, there is some arbitrariness in the definition 
of the massless fields. Concretely, we have adopted a basis in which the hidden photon does not 
couple to the SM particles]^ Hence the massless connector between the SM and hidden sector is 
the ordinary photon. This choice makes clear the fact that the dominant process to create hidden 
particles or to transfer energy from the visible sector to the hidden sector is through the creation of 
DM pairs, see Fig. [T] In this basis, by far the most important process to do so is DM creation from 



SM fermions through the photon, or Z boson, in the s-channel, // — )• e'e', see Eq.(72) of Appendix 
A. There is also the W^W~ — )• DM DM process but it turns out that it has always a small effect, 
even in regimes for which ^/s > 2mw- 

Note that, while in all the plots we take into account all contributions, in this section and in the 
next section our discussions will be focused on the case where the production of DM is dominated 
by non-resonant scattering processes. For instance, for // — )• e'e' processes, we put apart the 



^ The basis of Eq.(|3| is convenient for our purpose, as we focus on the possibility of creating DM through the 
kinetic portal. In particular we consider that the Universe contains initially only SM degrees of freedom, among which 
are ordinary visible photons. By definition the visible photon is the spin-one massless particle that couples, say, to 
the electron. In this basis, it is also clear that hidden radiation {i.e. extra degrees of freedom) may be produced only 
through the production of DM particles. We could have worked in a different basis though. In particular, instead of 
the transformation of Eq. ([3|, which gives 7 couplings to both e and e', and a 7' coupling to only e', we could have as 
well done a transformation which gives 7' couplings to both e and e', but a 7 coupling only to e. Naively we would 
then expect that the Compton scatterings ee — >■ 77' and 67 — >■ ej' play also a role, but this is not the case |13| . In 
any basis, the actual visible photon is by definition the state that couples to electrically charged SM particles. If this 
state is represented by a mixture, then the amplitudes must be added coherently. Hence, in this basis too, it is clear 
that our initial condition for the Universe contains only (essentially) visible radiation if there is initially no (resp. 
negligible) DM. Physically, it is also clear that the energy transferred to the hidden sector must involve the creation 
of a pair of DM particles, a process which is manifestly basis independent. 
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Figure 1: Processes that are relevant for the production of DM and thermahzation of the hidden 
sector through kinetic mixing. We work in the basis in which the DM is miUicharged so that the 
hidden photons, 7', only couple to DM and not the SM degrees of freedom (see text). 



contribution from Z — t- e'e' decay. These discussions therefore apply to the case of itldm > Tnz/2 
case, where obviously there is no Z — )• e'e' decay. It also applies to rriDM ^ 1 GeV. Indeed, in 
the latter case the Z decay contribution can be neglected because, as we will see, the production 
from the 7 mediated // — )• e'e! scattering is enhanced at low temperature (being maximum at 
T ~ Max[m£)Af , mj]), whereas the production from the Z decay is Boltzmann suppressed at these 
temperatures. There is however an intermediate mass range, 1 GeV < moM ^ for which the 

Z — )• e'e' decay may dominate the DM production. In this case, it is technically more convenient to 
express all Boltzmann equations directly in terms of the decay width rather than to hide it in the on- 
shell part of the scattering contributions. This we have done in section 6 for the Higgs portal because, 
in this case, the Higgs decay to a DM particle pair (if allowed, thus for rrioM < Tah/2) always 
dominates the production of DM. For a better comprehension of results within this intermediate 
mass range we therefore refer to section 6, where all equations apply in the same way, as well as 
to Appendix D, where a few supplementary subtleties concerning the interplay between decay and 
scattering in the case of the kinetic mixing portal are discussed in details. 

Note also that, a priori^ we would expect the scattering fe' — ?■ fe' of DM particles with SM 
fermions to be also relevant for the transfer of energy between the visible and the invisible sector. 
Ev However the behaviour of this t-channel process is very different from that of the s-channel 
processes. To begin with, the energy transfer in this t-channel process is a priori sub-dominant 
compared to that from the s-channel processes. Moreover this process involves a hidden sector 
particle in the initial state, whose number density with respect to the one of the visible sector 
particles is suppressed by a factor of with ^ = T' /T the hidden-to- visible sectors temperature 
ratio. However, the t-channel has a collinear divergent behaviour at low-energy transfers which 
could compensate for the two suppression effects. Nevertheless, we have checked that it is not the 
case and that the t-channel processes can be safely neglected in the calculation of the energy transfer 
between the visible and the hidden sector, as we will do in the sequel. 

Since the processes connecting the visible and hidden sector discussed above depend only on the 
combination 

K = e\/ a' I a , (5) 

in what follows, we will express all results in term of this parameter, k, which we call the connector 
parameter, and of a' , the equivalent of the fine structure constant in the hidden sector. The later 



coupling controls the process e'e' -H- 7'7' that takes place in the hidden sector, see Eq.(75) of 
Appendix A. 

To study the evolution of the number density of DM particles, we will solve a simple Boltzmann 
equation in presence of both the connector processes, SMiSMi o DM DM with i the various SM 
species, and the hidden sector process, e'e' -f-T- 7'7'. It is convenient to write this Boltzmann equation 
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in terms of the so-called 7 = Tueq reaction densities, defined in Appendix B. It takes the form 

dY / \ / Y'^ \ 

''^^ = E^™4i - +7H5(i - , (6) 

with s the entropy density, H the Hubble constant and Y the DM number density to entropy 
density ratio, Y = Ug,/ j s. Note that Ydm = 2y since DM is charged. Equivalently, this Boltzmann 
equation may be expressed in terms of the usual thermally averaged cross sections, (cu), with 
7 = {(yv)nl^ = Tneq, 

'J^ = T.^^connectv)^{Y^,{T) - Y^) + {a H Sv) {Y^,{T') - Y^) . (7) 

i 

Note that, in either form of the Boltzmann equation, we need to distinguish Yeq{T) = neq(T)/s from 
Yeq{T') = neq(T')/s. Thcsc abuudauccs involve the DM equilibrium number densities expressed as 
a function of T and T' respectively, while s is the total entropy density (visible plus hidden sectors). 
The first one, Yeq{T), parametrizes the number of SM particles participating in the SMiSMi — )• 
DM DM processes. The second one, Yeq{T'), parametrizes the number of 7' that are participating 
in the 7' 7' — )• e'e' processes. This distinction is important, as in the following we will consider 
situations for which T' <^ T, corresponding to Yeq{T') <C Yeq{T). 

If the DM particles never thermalize with the hidden photons, which, as we will see, actually 
happens if the combination Ka' is small, the Boltzmann equation for the DM abundance simplifies, 

^ / y"^ \ 

^""^^ = Yl ^connect (l " ^27^) " (8) 



■eq 



In other words, the number of DM particles is driven by the balance between DM pairs creation 
from SM particles and reciprocal processes, pretty much as in the standard WIMP framework. 

If, instead, the hidden sector thermalizes (in the sense that the hidden photons and the DM 
reach thermal equilibrium at some temperature T' < T) the {cfhsv) term must be taken into 
account. In this case one also needs to determine Yeq{T')^ and thus the evolution of the hidden 
sector temperature T' as a function of T. In other words, we have to determine the energy density 
of the hidden sector, /j', as a function of T, which is given by some possible initial hidden sector 
energy density, plus the contribution that may be transferred from the SM sector. The energy 
density transferred to the hidden sector by a process of the form 12 —t- 34 is given by the following 
Boltzmann equation 



dfj_ 
~dt 



+ 3Hip' + p') = I n^^^K •5i/i(Pi)/2(P2)KA^i2«34|'(2vr)V^)(pi + - - p^) ^Etr , (9) 
•' i=i 

which, as shown in Appendix C, may be rewritten as 

^ = -777^^ / ds . a{s){s - 4m')sTK,{4) , (10) 
where p is the energy density of the visible sector and K2 is the usual modified Bessel function. 



For T > m£,M {T < mDAi) a numerical resolution of Eq. (10) gives that p' / p increases like ~ 1/T 
(~ const respectively). This behaviour can be seen in Fig. [2j which gives p' / p as a function of the 
visible sector T, for various values of the connector coupling and DM mass m^Af 

^Notice that the wiggles for T > 10 GeV seen in Fig.|2]in the panels corresponding to the three lightest candidates 
are associated to resonant energy transfer through the Z resonance. For the two lightest candidates, this contribution 
has little impact on p'/P (note the log-scale). For the moM = 10 GeV candidate, the energy transfer is on the contrary 
dominated, and enhanced, by Z decay. See Appendix D and also section 6. 
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Figure 2: Evolution of the ratio of the visible (p) and hidden (p') sectors energy densities, for a 
range of connector parameter, k = 10~^'~^'~^'~^'~^'''~^^ (from up to down), and for various DM 
masses. 

Notice that for these calculations we have assumed that the DM is in thermal equilibrium with 
the 7' as, for all practical purposes, energy transfer is relevant only in this case. Indeed, if the hidden 
sector does not thermalize, the knowledge of p' is unnecessary in Eq. (j6| since both the Y^^{T') 
source term and the small contribution of the hidden sector to the Universe expansion rate may 
be neglected. Now the condition for thermal equilibrium of the visible and hidden sector through 
the connector is simply {crconnectv)neq ^ H which, taken at T ~ thdm, translates into roughly 
K > 10"'^'"'''"'''"^ for rriDM = rn^, 0.1 GeV, 10 GeV, ITeV respectively. For larger k, T' reaches T 
and p' /p stops increasing. As can be seen in Fig. |2]it may even decrease because, after the visible 
and hidden sectors have decoupled, reheating occurs in the visible sector each time a SM particle 
species becomes non-relativistic {i.e. T' /T decreases). We also notice that, below the critical value 
of the connector coupling at which the hidden sector may thermalize with the visible sector, p' /p 
reaches a plateau at T ~ moM^ which is given by p' / p ~ K?a'^mpi/{m£)M{g*^^ {T = rriDM))^^'^)- 
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3 Dark matter relic density: the Mesa phase diagram 



In the following we will first assume that the initial hidden sector energy density is negligible with 
respect to the energy which has been transferred from the visible sector at the temperature where 
the DM relic density freezes. At the end of this section we will discuss what changes when this is 
not the case. 

Assuming a negligible initial hidden sector energy density, the relic density is determined by 
three parameters only: the DM mass muM, the connector parameter k, and the hidden sector 
coupling a'. By construction there is no dependence on any other interactions or new particle 
masses, since in the visible sector we assume nothing but the SM. 

Fig. 3 shows the results we get for the relic abundance as a function of k and a' and for various 
values of the DM mass. The figure has the simple characteristic shape of an isolated flat-topped 
hill, or what we call here the "Mesa". The figures reveal the existence of essentially four distinct 
regimes or phases, which we will explain in details in the rest of this section]^ 

The first regime (phase I) is that of freeze-in, which corresponds to no thermalization, either 
through the connector or the hidden gauge interaction. Alternatively, if thermalization takes place, 
freeze-out of the hidden gauge (phase III) or connector (phase IV) interactions can occur, depending 
on which interaction is dominant. A fourth possibility, which we will explain in details, is due to 
the existence of an intermediate reannihilation regime (phase II). In this regime, there is a subtle 
inter-play between the connector and hidden gauge interactions. Like in the other regimes (freeze- 
in and out) we show that, in the reannihilation regime, the relic abundance is given by a simple 
analytical expression. 



3.1 Phase I: the freeze-in regime 

If both K and a' are sufficiently small, none of the interactions may thermalize the DM particle, 
neither with the SM sector, nor with the 7'. The relic density is therefore given by the (by now 
standard) freeze-in mechanism [2l |3]: the number of DM particles is simply given by twice the 
number of pair creations. This leads to the left-hand-side cliff of the Mesa shown in Fig. [3] the 
relic density is independent of a' since these interactions are negligible, and in this regime the relic 
density only depends on k. 

In practice it means that the last three terms of Eq. ([7| can be neglected and the relic density is 
simply given by integrating ^connect over time, which gives the number of e' particles produced per 
unit time per unit volume. Equivalently it is given by the integral over temperature of dY/dT = 
—7 connect / (T H (T) s). The freeze-in production is infrared dominated because, for large T, one has 
dY/dT ~ (for a cross section which behaves like 1/s for large values of s, as in Eq. (72)). 



As a result we can approximate the total number of e' particles produced by the number of e' 
produced per unit time, jconnect, times the Hubble time at T = max[mi,mDM], with rrii the mass 
of the initial SM particle in the pair production process. At later times the production becomes 
Boltzmann suppressed because less and less SM particles may produce DM pairs. In other words, 
Y grows monotonously as l/T until it reaches a plateau where it freezes-in, with the value 



-t r ^connect 

Y = c 



sH 



(11) 

T=max [»Tti,m£)jvf] 



*Note that in Ref. [4] phase diagrams have been considered for the more general case where there is an additional 
feebly coupled A particle in the visible sector that decays or annihilates slowly to DM particles. In this case, the 
phase diagram depends, in addition to the DM mass and the connector and hidden sector interactions, also on the 
mass of the particle A and its annihilation rate in the visible sector. The phase diagrams we have obtained, and 
shown in Fig. |3j are simpler and more predictive, as the visible sector, on which they are based, depends only on 
known interactions and parameters. 
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Figure 3: Phase diagrams for the kinetic mixing portal: contours of Y^m as a function of k and 
a' for rriDM = ""^e; 0.1 GeV, lOGeV, 1 TeV. The dashed thick hne gives n^Mh'^ = 0.11, or in 
other words YoMfnoM = 4.09 • 10"^'^ GeV. We have drawn the "transition hues" dehmiting the 4 
phases. Phases I, II, III, IV correspond to the freeze-in, reannihilation, hidden sector freeze-out and 
connector freeze-out phases respectively. There are transition lines between I and II, II and III, III 
and IV and I and IV. The solid black line corresponds to e = 1/ \/47r. Below this line the connector 
interaction is not expected to be perturbative. Note that, for rriDM = lOGeV, the Z boson leads 
to a "Mesa" phase diagram that is slightly more complex, but analogous to the one obtained in the 
case of the Higgs portal (see Section 6 and Appendix D). 

with c is a coefficient of order 1 (and with, as said above, Ydm = 2y). Numerically we get 
Ca = 5.0(4.8) for the e^e~ — )• e'e' channel and rriDM < fne (resp. ttidm > fne)- For — ?• e'e', 
we have c = 9.1(4.8) for uidm < "i^ (resp. ttidm > Trifj,)- The c coefficient depends on the s 
dependence of the cross section considered. For » "ij the actual value of rrii is irrelevant 

and the various fermions distinguish themselves only through their respective electric charge (and 
number of colors). The c coefficients are slightly larger than unity because the production occurs 
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Freeze-in (0.0^^1^=0 ■ I 




10-6 iQ-4 0.01 1 100 10" 

'«DM (GeV) 

Figure 4: Values of k that give the observed rehc density through freeze-in (a' is assumed to be 
neghgible). The continuous hne corresponds to the contribution of both the 7 and Z channels, the 
dashed line is only for the 7. 

mostly at a temperature that is a factor of 2-3 times smaller than max[mi,mDA4]- This is the 
temperature at which Boltzmann suppression is actually effective. 

Provided that the production of DM is non-resonant, analytically we get that the abundance 
of DM particles scales as F ~ 'yconnect/isH)\T=max[mf,mDM] °^ '^'^I'^DM if ruf < TUDM, and as 
Y oc K?/mf if mj > ttidm- The channels that dominate the production are those for which 
THf < rriDM, since they are effective down to T ~ ttidm- Therefore, if the candidate is heavier 
than the electron, Y is independent of m/ and scales as K^/niDMi and so 0_dm is independent of 
rriDM (modulo threshold effects each time mo hi becomes larger than the mass of one of the SM 
particles). This can be seen in Fig. |4j which gives the value of k required to have the observed DM 
relic density. If instead moM < me, then Udm scales like k^ijidm, as can also be seen in Fig. |4j 
Finally, for muM ^ 1 GeV and up to mz/^, DM production is dominated by a resonance, i.e. Z 
decay. In this case, the DM abundance is independent of moM (see Appendix D and also section 6) 
and scales like k^, so that Qdm scales like k^itidm- This behaviour corresponds to the dip in Fig.jij 
Notice that it implies that a smaller value of k is required to reach the observed DM abundance, a 
feature also visible in Figj3] 

The various scaling properties discussed above can also be seen in Fig. [5| which displays Y 
as a function of k. for various values of a' and rriDM- In this figure one clearly sees, for tiny 
values of a', the characteristic volcano shape of the transition between the freeze-in and freeze- 
out regime, from the processes driven by k. The top of the volcano corresponds to the point 
where the connector interaction thermalizes, delimiting the ^dm ~ freeze-in behaviour from 
the 0,£)M ~ ^/ {c^connectv) ~ freeze-out behaviour. In Ref. j3] similar transitions have been 

obtained for other types of interactions. For larger values of a' however, there is no more freeze-in- 
to-freeze-out transition. The volcano becomes a truncated volcano and the transition from freeze-in 
to freeze-out undergoes an intermediate regime of reannihilation, which we will now discuss. 

3.2 Phase II: freeze-out with a source term: the reannihilation regime 

Starting from a freeze-in situation, if one increases k and/or a', at some point, to be defined below, 
there are enough DM particles and the interactions in the hidden sector are fast enough for the 
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Figure 5: DM relic abundance Y^m as a function of the connector parameter k for different DM 
masses uidm and values of the hidden sector interaction, log^Q(a'/a) = 1, —1, —3, —5, —7 (bottom- 
up). 



hidden photon 7' and the dark matter particle e' to thermalize, even if the connector interaction 
remains out-of-equilibrium. Therefore one may in principle define a hidden sector temperature T' , 
with T' < T. In this case, we may determine the dark matter relic density in two steps. First we 
estimate the hidden sector energy density, p', as a function of the visible sector temperature T, as 
explained above, and define a hidden temperature T' , which depends only on k and moM, through 
p' oc T'^. This step gives us in turn a way to define the DM equilibrium number density in the 
hidden sector, Yeq{T') 7^ Yeq{T). From these, in a second step, we may integrate the Boltzmann 
equation for the DM number density, Eq. ([g]), taking into account both the source connector and 
annihilation in the hidden sector. 

Starting from an initially insignificant abundance of hidden sector particles, the DM number 
density follows the freeze-in regime until the 7' and e' thermalize. The condition for chemical 
equilibrium between the 7' and e' is 

^annih = {cT H Sv)neq{T') > H , (12) 

whereas kinetic equilibrium holds provided 

1 dp'/dt 



H p> 



>1, (13) 



where p' here stands for the energy transferred from the DM particles to the 7'. Actually both 
conditions essentially coincide because the DM and 7' thermalize for T > rn£)M when the DM 
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particles are still relativistic. In terms of the cross sections, Eq. (12) can be written as 



{aeffv)neg{T)^/c> H with {(Jeffv) = \J {(J HSV) {(J connect) , (14) 



where c is a constant of order unity defined in Eq. ( 11 ). This reflects the fact that neq{crvHs) is pro- 



portional to the number density of DM particle, which itself is proportional to {aconnectv)n'^g{T) / H . 
Taking this condition at T ~ rrioM gives the "phase transition" line between the freeze-in and 
reannihilation regimes. It gives a value of a' larger than a quantity which scales as oc moM / {itipik) 
(numerically we get Ka' > 1 • 10~^^GeV~^ ■ fnoM /{g*^^ {T = rnoM))^^^ with g^^-^ the effective num- 
ber of degrees of freedom of both SM and hidden sectors). Once the hidden sector has thermalized, 
the DM number density is just that of equilibrium, which for T' > ttidm is 

Y.,(r) = '-^^f. (15) 



g*s 



in which ^ = T' /T < 1, ge' = 2 is the number of degrees of freedom of e' and T' is defined by the 
equilibrium relation 

P='^9*'T'\ (16) 

Later on, the number density Y follows the equilibrium one Yeq(T') and does so until the latter 
becomes sufficiently Boltzmann suppressed. However the Boltzmann equation at this stage is differ- 
ent from the one in the standard freeze-out as it contains an extra source term from the connector. 
This makes the reannihilation regime a bit complex, so let us explain what happens step by step. 

The reannihilation regime corresponds to a situation in which interactions within the hidden 
sector are fast (that is, compared to the expansion rate), while the energy transfer from the visible 
sector is comparatively slow. In this case, the hidden sector reaches thermal equilibrium, but at a 
lower temperature than that of the visible sector, ^ = T'/T <C 1. Hence the particles in the visible 
sector are more abundant (per species) than in the hidden sector, Y <C Yeq{T). Starting from the 
full Boltzmann equation Eq. ([T]), we thus have 



HdY 

s dz 



{aconnectv){Y,liT) - Y^) + {o H sv) {Y!^{T' ) - Y^) 

{cTconnectv)Y^q{T) + {a H Sv) {Y!q{T' ) - Y^) , (17) 



(here we have dropped the initial SM particle i index to avoid the cluttering of symbols). 

By assumption, early on the processes DM DM o j'^' are fast, and the abundance of DM is 
able to track Yeq{T'), but when T' < rriuM the equilibrium abundance Yeq(T') becomes Boltzmann 
suppressed and so does the rate of ■j''y' — )• DMDM. If the connector source term is negligible, 
decoupling takes place when Yeq(T') drops below the critical value 

Y,rit = H/{CJHSV)S ~ Y,q{T',„^) , (18) 

a condition that corresponds to 

{oHsvWqiT') ~ H{T). (19) 

This is just like in the standard freeze-out mechanism. 

If, however, the connector source term is large enough, then, at some temperature > Tcrit, 
the 7' 7' — )• DMDM process may become sub-dominant compared to the production of DM through 
the connector. This occurs provided 

{<y^nnectv)Y^q{T) > {a H Sv)Y^q{T' ) (20) 
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for T <Ta- Clearly, the connector source term may become more important than the hidden sector 
source term simply because T' < T, so that Yeq{T') <^ Yeq{T) for T' < moM- This is illustrated in 
Fig. |6] (in particular the right panel), where l^g(T') becomes Boltzmann suppressed while Y(,q{T) ~ 
const {i.e. T > ttidm > T'). In this situation the 'y''j' — )• DMDM process becomes irrelevant 
around Ta, and this is true even if its rate is still larger than the Hubble rate. Consequently, from 



Eq. (17) we get that, for T < Ta, the DM number density is determined by the following simpler 
Boltzmann equatioDj^ 



dY {(Tconnectv)s 2 (rp-, _ {'^HSv)s 2 

dz ^ H '^^^^ H 

{(^ connect^) S ^^2 frr\ ^ 

crit 



rr Ye^T)-^. (21) 

±1 1 rt 



Now, down to T ~ T^, Y ~ Yf>q{T') holds, as the rate for DM pair creation by hidden photons 
is faster than the Hubble rate (and the rate of production through the connector). Moreover, at 



T ~ Ta, first Y > Ycrit holds, and second, using Eq.(20), one has 



(T) ^ ^^MgHy^ (22) 



H ' H Y, 



crit 



SO that both terms on the RHS of the Boltzmann equation (21) are relevant at T ~ ^aj^ Hence 
below Ta the tracking solution is given by y = Yqse, 



—V {'^ connectv) S 2 ( rp\ _ {'^ connect ^ / ^2 / rp\ /r,n\ 
^QSE = J^crit ^ J^egU ) — ~/ ^-'egU ) ' K'^'^) 



that is dubbed the Quasi Static Equilibrium (QSE) abundance 14] to emphasize its intermediate 
character, see Fig. [6j 

Yeq{T') < Yqse < Y,q{T) . (24) 
The abundance of DM tracks the QSE until a temperature T = Tf at which Yqse drops below 



the equilibrium condition critical value of Eq. (18). Indeed at T = Tj the condition Yqse — Ycrit is 



equivalent to the conditions (aconnectv) sY^q{T) / H ~ Y and Y ~ Ycrit, which shows that both terms 
decouple at the same temperature Tf. From this point on the abundance of DM is frozen. 

In other words, and to summarize, the condition that has to be satisfied to stay in the rean- 



nihilation regime is that Eq. (20) holds before the freeze-out of interactions in the hidden sector, 
Eq. (19). If it is not the case, we have a standard freeze-out scenario. Otherwise, as illustrated in 
Fig. |6| the abundance of DM begins to follow the equilibrium density of the hidden sector Yeq{T') 
until it intercepts at a temperature T = Ta the QSE abundance, which becomes the new track- 
ing solution. Equilibrium between the connector source term and annihilation into hidden photon 
freezes at T = Tj, and the final DM abundance is given by the value of the QSE abundance when 



it intercepts the critical density of Eq.(18), 



Y{Tf) = YcruiTf) = YQSEiTf) . (25) 



From the Boltzmann equation (23) we may get an analytical estimate of the final relic density. 



First using Eqs. (18) and (23), we get from Eq. (25) 



Y\TA = ^"""^'^' YliTf) . (26) 

{(JHSV) 



similar equation has been considered in Ref. [4] for the case of a decay process (rather than annihilation) as a 
source term from the visible sector. 

*That is to say, each term gives in the Boltzmann equation a contribution to {z /Y){dY / dz) that is larger than 1. 



13 



Region.II mDM= 200GeV, /<:=10"*, ff' = 10"^*) Region.II (m^f^^ 200GeV, /f=10"'^^ (i' = 10"^'') 




0.1 I 10 100 0.1 1 10 100 



z='«dm/T z=mDM/T 

Figure 6: Examples of evolution of the DM number density y as a function of z = niDM/T, in 
the reannihilation phase (red hue). Also shown are Yqse, Ycrit, Yeq{T') and Yeq(T). The left panel 
is obtained for tudm = 200 GeV, k = 10"^ and a' = lO'^-"^. It gives T'/T ~ 0.067 at T = 
and T'/T ~ 0.072 at T = Tj. The right panel is obtained with rriDM = 200 GeV, k = 10~^'^ and 
a = cases lead to a relic density that is in agreement with observations. 



that, using Eq. (18) once more in Eq. (26), gives 



neq{Tf)^/ {(Tconnectv){(JHSv) = H{Tf) 



(27) 



This condition is like the one for standard freeze-out, but with an effective cross section given by 
{(Teffv) = y/ {(Jconnectv) {cfhsv) Solviiig for Tf in the usual way we get 



9e 



log[0.038^^=mpimDM {cTef fv)c{c + 2)] 
1 



9e 



log[log[0.038 

2 Igeff 



mpimDM{creffv)c{c + 2)]] 



and 



YeqiTf) 



3.79x fygl 
gi,smpimDM{creffv) ' 



(29) 



(30) 



with Xf = mohi/Tf and c a numerical constant of order unity. Using Eq. (26) we get finally 



Y{Tf) = Yqse{Ts) 



3.79xfy 
g*smpimDM{crHsv) 



^Equivalently, for Y ~ Yqse we may rewrite the Boltzmann equation ( |2l[ ) as 

H dV 



Z — ~ 2 a/ {Oconnectv){oHSV)Yeq {YqsE — Y) , 

s dz 



(31) 



(28) 



which shows that freeze-out takes place at a temperature T — Tj such that H ~ \/ {o'connectv){o'HSi^)neq{T) and 
Y = YQSEiTf). 
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For the standard choice of numerical constant c(c + 2) = 1, this expression approximates the exact 
result within a factor less than 2-3 in the reannihilation regime of Fig. [3] Better accuracy may be 
obtained by adopting sligthly smaller values, for instance c = 0.1 for moM = 1 TeV. Note that Y 
depends on {(Tconnectv) only through the logarithmic dependence of xj. 

To sum up the final relic density obtained in the reannihilation scenario is inversely proportional 



to the hidden sector cross section, Eq. (31). This is on one hand the standard expectation. On 
the other hand its decoupling temperature, x/, is determined by an effective cross section (aeffv) 
which is the geometric mean of both connector and hidden sector cross sections. In practice since 
the relic density depends linearly on Xf, this can change the result by up to one order of magnitude. 
For example with a hidden sector cross section of the order of the one that is typically needed in 
the standard freeze-out mechanism, xj can be as low as ~ 2-3 instead of the usual ~ 20. In this 
circumstance, to get the observed relic density, the hidden sector cross section has to be one order of 
magnitude smaller than in the usual freeze-out scenario. Note also that xj cannot be smaller than 
~ 2-3 because the source term, and therefore Yqse, gets Boltzmaim suppressed when T < uidm- 
Also it cannot be larger than the usual ~ 20 value because, using the observed relic density in 



Eq. ( |31[ ), one gets anyway a value of (ansv) which, plugged in Eq. (29), cannot give a larger value 
of X/ (as a connect < ctrs in this regime). As for x'j = ttidm /Tj it cannot be smaller than ~ 2- 3 
too because Yeq{T') decouples when T' < itidm, but it could be orders of magnitude larger than 
unity (if T'/T" << 1)- In all cases DM freezes when non-relativistic because the DM particles left 
mainly consist of particles which have been pair created at T ~ Ty with Tf < m dm (the particles 
produced early on having annihilated to 7' since a while) and because anyway, when they freeze, 
they are still in kinetic equilibrium with the 7' and this down to a temperature T' which is below 
T' 

From the results above, the structure of the phase diagram in the reannihilation regime can be 
easily understood. For a fixed value of the connector k, as a' increases, il-DM obviously decreases 
because, as usual, Y{Tf) oc l/{aHsv) decreases. For a fixed value of a' the dependence on k on 



the other hand is milder because it enters only logarithmically, through xj, Eq. (29). In the phase 
diagram, this explains the near flatness of the part of the Mesa corresponding to the regime II in 
Fig. [3] 

We note also that, in practice, if the DM production is dominated by processes where the 
mediator mass is irrelevant (as applies obviously for ttidm > ^zf^, or for moM ^ 1 GeV where 
the 7 mediated processes dominate), we did not find any case for which, without thermalization of 
the connector, the DM particle would freeze-out without experiencing a period of reannihilation, 
i.e. where Y follows the QSE distribution for a period of time. This can be understood in the 
following way. Starting from y = 0, if the hidden sector thermalizes, i.e. if Y goes over Ycrit, Y 
increases quickly up to a point at which Y = Y(,q(T') > Yqse > Ycnt^ Later on this hierarchy is 
not modified down to T' < m/jA/ because the three abundances have similar dependence in At 
T' ^ rriDM, Y becomes rapidly Boltzmann suppressed while Yqse (which is Boltzmann suppressed 
only at T < m£,M) and Ycrit do not change their behaviour. As a result Y intercepts Yqse before it 
intercepts Ycrit and reannihilation occurs before freeze-out. Only on the border line between freeze- 
in and reannihilation does DM just reach the critical value before immediately decoupling (without 
knowing either a sizable period of reannihilation or a sizable period of thermal equilibrium). 

The range of parameters for which the reannihilation regime actually occurs is quite large, see 



This can be anticipated from Eq. (231, wliicli at the thermalization point gives Yqg^ ~ YcritY , i.e. Yqse = Yc 
(since, before thermalization, Y ~ {a'connectv)sY^q{T)/ H), and the fact that just after thermalization {Y > Ycr 
Yqse still increases essentially as the geometric mean of Y and Ycrit- 



9, 



fHsu) is essentially constant during this period because T' never go much above ttidm- Therefore from Eq. ( 18 1 



we get Ycrit ~ 1/T and from Eq. (231 we get Yqse ~ l/T. As for Y^qiT'), using the p'/P ~ ^JT relation above, it 
scales as Ycq{T') = ncq{T')/s ~ T'TT^ ~ l/T^/". These scaling properties can be seen in Fig.je] 
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Fig. [3] In particular, values of parameters which could lead to a direct detection rate compatible 
with, for instance, the results of the DAMA or CoGeNT experiments (see Section 4) turn out to lie 
in the reannihilation phase. Note also that the period of reannihilation can cover several orders of 
magnitude in T. 

3.3 Phases III &; IV: Freeze-out with hidden gauge or connector interactions 

Phase III of the Mesa in Fig. [3] corresponds to the situation in which both the hidden sector and 
connector interactions are fast enough to thermalize but in such a way that the hidden sector 
interaction annihilation rate remains dominant in the Boltzmann equation. In other words both 
sectors share the same temperature, T = T' , but once they have thermalized, the connector not 
only does not lead to any further energy transfer between both sectors, but also play no role in 
the freeze-out of DM, which is thus determined only by a' . This explains the shape of the phase 
diagram for large values of a' (top of the Mesa, where the relic density is independent of k for fixed 
a'). In this case, the relic density simply scales as Qdm oc 1/ {(Jannih.v) oc mj^j^.j / a'"^ as can be seen 
in Fig. |3] and [5] Also, in this case there is no reannihilation processes since the source term can be 
neglected as soon as T ~ T'. 

The "phase transition" line between the reannihilation and hidden sector freeze-out regimes 
corresponds to the thermalization condition {(7connectv)neq{T) / H\T~mDM > considered at the 
end of Section 2. This condition applies also for the freeze-in to connector freeze-out regime phase 
transition, relevant for very small values of a'. Numerically we get k^8 ■ 10~^ ■ ruj^'^j / {gt^ ^ {T = 

Finally, if the connector parameter is so large that it does all the job of thermalizing both sectors 
and also dominates the DM freeze-out, then, obviously, the hidden sector interaction becomes 
irrelevant, which explains the vertical behaviour of the phase IV of the Mesa in Fig. [3] with a 
relic density which decreases as the connector increases, as usual through standard freeze-out, 
^DM ^/ {(^annih.v) OC m?j^j^^/K? (as Can also be seen in Figs. [3] and [s]). In this case a' plays 
little role, apart from creating a thermal population of 7'. The phase transition line between both 
freeze-out regimes simply stems from the condition {aconnectv) > {(Jhsv)- Numerically we get 
K > 4.6 • 102 . a'/ (<7f^^(r = mDM) f^. 

Note that, for fixed values of a' and increasing values of k, the reannihilation phase goes nec- 
essarily towards the hidden sector interaction freeze-out phase, rather than directly towards the 
connector freeze-out phase. This stems from the fact that in the reannihilation phase the hidden 
sector interaction is already in thermal equilibrium whereas the connector one is not. The border 
between the reannihilation and hidden sector freeze-out phase corresponds to the stage where the 
connector is just large enough to begin to thermalize, whereas the hidden sector is already in deep 
equilibrium. Therefore in this case the connector decouples before the hidden sector interaction 
does. 

Note also that all phase transition lines meet at a single point. This is due to the fact that at 
the point where the II-III and I- IV line (where T connect = H) meets the I-II line (where Ths = H), 
one has necessarily Tconnect = ^HS which is the condition the III-IV line fulfills. At the meeting 
point both cross sections decouple as soon as they thermalize. 

3.4 What if the primordial hidden sector energy density does not vanish? 

As can be seen from the above discussion, to assume an initial population of DM particles in 
the hidden sector, and to neglect the effect of any possible connector {i.e. considering T' /T as a 
constant), as has been done in previous works (see Refs. |10 | [TT | [T2]). is justified only if the connector 
is very tiny, actually even smaller than in the pure freeze-in scenario. If this is the case one is left 
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Figure 7: Values of a' required to get the WMAP relic density as a function of moM assuming 
no connector between the hidden sector and the SM sector, for different values of the temperature 
ratio ^ = T IT = 0.01, 0.1, 1 (bottom-up). 

with an ordinary freeze-out scenario in the hidden sector, induced by the hidden sector annihilation, 
but still mediated by a Hubble expansion rate that receives contributions from the abundance of 
both the visible and hidden sectors particle species. The relic density depends on the parameters 
TTie', a' and ^ = T' /T. In the non-relativistic limit, the s-wave, dominant part of the cross section 
of Eq. ( 75 ) (reversing the in and out particles in this process) gives {(Tannihv) 



DM- 



For 



the freeze-out temperature Tf we obtain: 

Xf = ^ ■ ^n[0.038 • f^'^{(Tannihv)mpimDM p — c(c + 2)] 

-^^ • • ^n[0.038 • {(7annihV)mpimDM-J^=c{c + 2)]} 

2 leff 



9* 



which gives 



2 1.07 X lO^x/ • GeF" 
ilDMh = 2 

{9*s/^glf^^--'- • 



(32) 



(33) 



mpi{aannihv) 



where here too a factor 2 comes from the fact that DM is composed of both a particle and its 
antiparticle. c is a numerical constant of order unity. We get that the values c = 0.22, 0.35, 0.4 give 
the correct abundance for ^ = 0.01, 0.1, 1 respectively. 

Here we have assumed that the present value of the ratio of temperatures, or better is small, 
so that (today) ~ 3.91 (which implicitly means that the Hubble expansion rate is dominated 
by the visible sector particle species). In Fig. [7]we show the values of moM and a' which, from a 
numerical integration of the Boltzmann equation, lead to a relic density of iloMh'^ = 0.11. This plot 
can be understood in the following way. Given the Boltzmann suppression of the DM equilibrium 
number density around freeze-out, riDMiT') oc (T')^/'^e~™^*^/"^ , the value of rriDM/T'^ = Xf does 
not change much with ^, even though the smaller ^ is, the larger is the Hubble parameter at T' ~ Tj- 
(when Tannih — H): the dependence is only logarithmic. In other words Xf = S^x'j changes linearly 



in ^, up to a logarithmic factor, see Eq. (32 ). But if one considers a smaller ^ the value of the Hubble 
constant (dominated by the visible sector, i.e. H ^ const ■T'j) at freeze-out changes as since 
H{Tf ) oc r? oc T'? 1^ . This implies that the number of DM particles left at freeze-out, UDM^'t)-, 



increases as (since H{Tf) = TannihiTj-) ~ {(Jannihv) np)M{Tjr))- However the entropy at freeze- 
out (also dominated by the visible sector, s ~ const -T^) increases similarly as Therefore 
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^DM = mDMn^D^^ / pToday ^ mDMiriDMiT'jr) / s{T'j-)) .{sxoday / PToday) IS proportional to ^, in the 
same way as Xf, see Eq. (33). For smaller ^ this must be compensated by taking a smaller (Tannih-, 
as shown in Fig. [7j so that freeze-out occurs earlier. Note that these results agree with the one of 
Ref. |TT| up to the power 5/2 instead of 3/2 in Eq. ( [32^ (which is freeze-out prescription dependent 
and has a moderate numerical effect). 

The condition that determines which DM source dominates the relic density today, the pri- 
mordial one or the connector induced one, reduces essentially to a condition on energy densities. 
For instance, taking values of parameters leading to the observed relic density in Fig. |2] an initial 
hidden sector population has essentially no influence if the initial p' / p\init. is smaller than the value 
of p' / p\Tf which, starting from zero initial abundance, one obtains at the temperature where the 
abundance freezes (in or out). Conversely, if this condition is not satisfied, the connector can be es- 
sentially neglected for the determination of the relic density. Therefore to know which contribution 
dominates, one can essentially compare the initial p' / p with the value one obtains from Fig. [2] It 
is a priori unlikely that both contributions would be of same magnitude. 



4 Testing freeze-in with direct detection experiments? 

An interaction between dark matter and quarks mediated by a massless gauge boson is interesting 
for DM direct detection in many respects. To begin with, it allows to give an observable signal 
even for very tiny value of the coupling because the cross section has a collinear infrared divergence 
leading to a direct detection elastic cross section proportional to i-e. which is enhanced at the 

low nuclear recoil energies Er- More mundane dark matter particles, which are supposed to interact 
with a nucleon through a massive particle (or at least heavier than the recoil energies which are 
typically in the few keV range, see e.g. |15| in the case of massive Z'), have a constant cross section 
on quarks. This is also a feature which is absent in other feebly interacting scenarios. A corollary 
of the energy dependence of the elastic cross section is that the recoil energy spectrum is distinct 
from that of generic DM candidates. In particular, this behaviour favours experiments with a low 
recoil energy threshold, a feature which has been invoked as a way of reconciling for instance the 
DAMA/LIBRA measurements with the exclusion limits by other experiments in the framework of 
the mirror model for dark matter |16| I17j . Notice that the velocity distribution of dark matter in 
the context of mirror dark matter is distinct from that usually assumed for WIMPs. Here we are 
conservative and make the standard choice of a Gaussian distribution with mean velocity that we 
set to vq = 220 km/s at the position of the Sun, truncated at an escape velocity Vesc = 550 km/s 
(see for instance jl8j). 

In this section we consider a few benchmark experimental results and confront the model with 
the data. The main conclusion is that present and future experiments may probe an interesting 
part of the parameter space of the model. In particular we show that Xenon IT will probe all the 
regimes of the model, including the freeze-in, a feature that is specific to the present scenario. Also, 
as already alluded to in the previous section, the reannihilation regime corresponds to candidates 
that may explain the DAMA [19j or CoGeNT data (see |20| I21| and also the discussion below), but 
these solutions are already strongly disfavoured by the current exclusion limit set by XenonlOO |22) . 
The main results are summarized in Figs. [8]and|9j 

The spin-independent (SI) elastic cross section between dark matter and a nucleus of mass 
number A and atomic number Z is, for kinetic mixing, given by 

-FiiqrA) , (34) 



dEr ^2^2 



■niA 
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Figure 8: Regions in the plane log]^o('*) ~^DM compatible with DAMA modulation data (red/dark 
grey blob) and CoGeNT total rate (orange/light grey finger) at 99 % C.L., together with the 
exclusion limit at 90 % C.L. from XenonlOO (see text for details). 



with Fa the Helm nucleus form factor 



qrA 



(35) 



q = (2mA-E'r)^^^ is the momentum transferred, va = 1.14^^/'* fm, s = 0.9 fm and ji is the n = 1 
spherical Bessel function of the first kind |17| . The cross section da/dEr may be cast in a form 
which is reminiscent of the SI cross section used by direct detection experiments 



dcT 
dEr 



niA 



[f^Z + UA-Z)f 



Jp 



with 



2 f -2 

f^nPM Ottk^Q;^ ■'P 

m\Er [fpZ + UA 



(36) 



(37) 



and where finDM is the nucleon/DM reduced mass and are the coupling strengths of p,n to 
the mediator. However this is of little practical use because, unlike in most models with a heavy 
mediator, Cgjj depends on Er- Hence we rather present the results in the k — uidm plane, with a' 
fixed by the requirement Q.DMh'^ = 0.11. 

These are given in Figs. [8] and [9] In the first figure we show, for the sake of illustration, a zoom in 



the region of parameters compatible with the DAMA/LIBRA modulation at 99 % C.L, The best 



fit is for K ~ 2 • 10 ^ and uidm ~ 15 GeV. A glance at Fig. [s] reveals that the candidates consistent 

^''To determine the DAMA/LIBRA region, we have used the detector resolution given in Ref. [23] . Also we have 
assumed no channelling, and the most standard choice of quenching parameters (concretely we have taken qi = 0.09 
an gjva ~ 0.3). The contour is given by a with n = 10 degrees of freedom (best fit point corresponding to 
K = 1.7 ■ 10"** and moM = 16 GeV, with = 5.2). 
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Figure 9: Exclusion limits at 90 % C.L. from the current XenonlOO data and forecast for XenonlT 
for one year (dashed, orange) and 4 years (dashed, green) exposures. The light grey line at the 
bottom corresponds to the pure freeze-in regime. The dip is due to the effect of the Z resonance 
on freeze-in. The upper, light grey lines delimit the region where reannihilation (below it) and 
hidden sector freeze-out (above it) regimes are possible (imposing the DM relic density constraint). 
There is a further division in the regime of hidden sector freeze-out (III), depending on whether the 
connector thermalizes (IIIB) or not (IIIA), see Appendix D. 

with both DAMA and WMAP are in the reannihilation regime II (actually close to the freeze-out 
regime IIIA). As Fig. [s] shows, they are also excluded by Xenon 1 00 For further illustration, we 
also show candidates corresponding to a recent, albeit preliminary, reassessment of the CoGeNT 
unmodulated data, as apparently a dominant part of the signal may at the end be attributed to 
As our purpose here is not to explain/reconcile the current direct detection 
experiments, we do not refer to the recent CRESST and modulated CoGeNT data. 

The bottom-line of this section is that currently, for rxiDu > few GeV, both freeze-out regimes 
are excluded, whereas the reannihilation regime is strongly constrained by XenonlOO. The freeze- 
in regime is allowed for any DM mass. As for the future, we refer to the Xenon-IT experiment. 
Assuming a one ton-one year exposure, and no event, the orange dashed curve of Fig. |9] gives the 

^^To derive our exclusion limit we have used a simple fit to the mean scintillation efficiency extrapolated to low recoil 
energies, as given in Fig. 1 of [22]. Concretely we have set Leff = for Enr < 1 keV, L^ff = —0.154 + 0.221ni5„r 
for 2 keV < E„r < 3 keV and Le// = 0.053 + 0.0353 In for Enr > 3 keV but below 100 keV. As in [22] we have 
taken into account the Poisson fluctuation in the number of photo-electrons (PE) events and have set the threshold 
at 3 PE. To set the limit we have used the Monte-Carlo code for the Optimal Gap Method, provided by Yellin and 
described in details in [24'. 

^^See the talk of CoGeNT at Taup2011. A substantial fraction of the original (unmodulated) CoGeNT signal 
[201 121| is now being attributed to background, so that the signal is much reduced, and also less peaked at small 
recoil energies. The residual signal is no more excluded by the CDMS low-energy analysis |26| . but it also moves 
away from the preferred DAMA region, making it more difficult to reconcile both experimental results, which is 
anyway not our purpose in the present work. Moreover, all the CoGeNT (tentatively reevaluated data) are excluded 
by XenonlOO. 



background [25] 
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sensitivity reach of XenonlT in the k — ttidm plant ^'^ Interestingly, this experiment may probe a 
large fraction of the reannihilation regime. Furthermore, it may also be possible to test the freeze- 
in regime. For the sake of illustration, we draw the exclusion limit that would correspond to an 
exposure of 1 year (orange dashed line) and 4 years respectively (green dashed line). In the first 
case, it could be possible to probe the freeze-in mechanism for a DM mass range between about 45 
GeV and about 100 GeV. With 4 years, the range extend to a candidate with a mass of about 500 
GeV. 

Positive signals with a recoil energy spectrum that is in accordance with a l/E"^ scattering cross 
section on nuclei would allow to distinguish this model from a more standard DM candidate, which 
usually predicts a constant cross-section. In principle, it should be possible to distinguish this model 
from mirror models, which a priori display the same l/E"^ dependence (see for instance |38)). but 
which have a distinct velocity distribution and, in general, a multi-component halo of dark matter 
with particles masses in the few GeV range. 

5 Compatibility with cosmological constraints 

Various cosmological constraints on DM particles interacting with massless gauge bosons have been 
considered in the literature |10| W\\ [121 123 [Ml [Ml ISO] . In this section we review the constraints that 
are the most relevant for our purpose, referring to the original literature for more detailed discussions 
if necessary, and extend some of the constraints to the case of a massive hidden photon]^ We show 
in particular how the stringent galactic ellipticity constraint gets considerably relaxed when one 
considers a slightly massive, rather than a massless, 7'. 

Primordial nucleosynthesis constraints: To begin with, let us dispose right away of the most basic 
constraint, based on Big Bang nucleosynthesis (BBN) and the increase of the expansion rate of 
the Universe in presence of additional relativistic degrees of freedom around T ~ 1 MeV (see for 
instance |29)). This is usually phrased in terms of limits on the number of neutrino families j31| [32] 
which, at 95 % C.L., we may take to be AA^^ < 1.4. In our simple model, there are only two extra 
degrees of freedom around the time of BBN (expect for the marginal case -i.e. in our work- of a 
light DM candidate, tudm ^ i^ie), corresponding to the polarization states of the hidden photon, 
so that 

£,^ = 0.16AN^ (38) 



p 10.75 

where, we recall, ^ = T' /T. BBN nucleosynthesis gives then ^ < 1.05. If the DM is lighter than 1 
MeV, the bound is a bit stronger ^ < 0.8 but still irrelevant for our purpose. A CMB constraint 
on the number of neutrino families Ni, < 4 leads to a similar conclusion |28] . The Planck satellite 
experiment is expected to narrow down the range to ANi, < 0.2. 

Galactic dynamics constraints: As one may expect, the most important constraints on our scenario 
rest on the (milli-) charged character of the DM and the associated long range interaction. In 
particular, this interaction may affect the way dark matter clusters in the universe, with potential 
consequences for the formation of clusters of galaxies, including the dynamics of the so-called Bullet 
Cluster, and modifications of the dark matter halo profile in galaxies. It turns out that the latter 



'^^For XenonlT, we have applied the same prescriptions as for XenonlOO regarding scintillation efficiency, the 
Poisson statistics for low PE and a threshold at 3 PE. 

^"^Notice that if we break the (7(1)' symmetry (either spontaneously or through the Stueckelberg mechanism), we 
lift the degeneracy between the visible and the hidden photons. In this section, we work directly in the basis of mass 
eigenstates. 
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effect gives tlie strongest constraint on tlie interactions between DM particles, and so on ol 
A long range force may affect the shape of dark matter halo in galaxies and clusters of galaxies. In 
particular they may erase the observed/inferred tri-axial elliptic shape of halo of elliptic galaxies. 
They may also lead to the formation of a cored (a central region with almost constant DM density) 
instead of a cuspy profile. The ellipticity of elliptic galaxies gives an upper bound on oi which 
is quite stringent [11]. In this reference, the bound is estimated to be given approximately by 



c^' i$ 10^"^ (m£)Af/GeV)^/^ (modulo an extra logarithmic dependence on rriDM and a'). Comparing 
this constraint with the values of a' reported in Figs, [sjand 10 (see also Fig. 11 ), we observe that the 
freeze-in regime, as well as the connector freeze-out regime is always allowed, since they work with 
any values of ol smaller than a given value which depends on the DM mass. But the reannihilation 
and freeze-out (see Figs. 3, 7 and Fig.l of Ref. [IT]) in the hidden sector regimes both require a 
fairly massive DM candidate, estimated to be rriDM ^ 1 TeV and rriDM ^ 10 TeV respectively. We 
notice however that the ellipticity constraint derived in |11| is tentative and that more work could 
change the limits on a' and ttidm- For instance, if we relax the ellipticity bound on a' by a factor 
of 3, reannihilation and freeze-out in the hidden sector are viable regimes for nioM ^ 100 GeV and 
nT-DM ^ 1 TeV respectively. 



LogijQ-'] (nDMh^=0.11) 
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Figure 10: Upper limit on k obtained from the ellipticity bound (solid line) on a' imposing that 
^DAlh'^ = 0.11, as a function of tudm- The dashed line gives the bound we get if we relax the 
ellipticity bound on a' by a factor 3. The lower line gives the value of k required along the freeze-in 
regime. This regime, ellipticity constraint included, allows any value of moM since it works with 
any arbitrary small value of a' . The contour lines give the value of a' one gets imposing the relic 
density constraint, in accordance with Fig. |3] 

Similarly the ellipticity of DM halos may put bounds directly on the coupling between DM and 
ordinary matter, i.e. on k. From Fig. 1 in |27| . we see that the constraint is mild, k < 10^^ for 

'^^The CMB puts very mild constraints on our model. For all practical purposes our candidates behave as cold dark 
matter and are essentially decoupled at the time of matter-radiation equality (this is not the same for the so-called 
Mirror Dark Matter, see for instance |29) . 
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iTT-DM ~ 10 GeV. Other effects [27], in particular tlie requirement of decoupling at recombination 
are more stringent, k < 10~^ again around uidm ~ 10 GeV. However, these constraints are weaker 
than the exclusion limits set by direct detection experiments ( |27] for CDMS-Si and Fig. [9] of the 
present manuscript for XenonlOO). 

Possible depletion of DM in the Galaxy: Direct detection constraints rest on the assumption that 
the energy density of DM in the vicinity of the Sun is pDM ~ 0.3 GeV/cm~^. This might not be 
the case for milli-charged particles, as has been emphasized in [33j and further discussed in |27j (see 
also [23])) as magnetic fields in the Galaxy may lead to a depletion of the DM abundance at the 
position of the Sun jS^. To be sure, this effect does not per se exclude a candidate. Rather it puts 
limits on its potential relevance for DM search through direct and, implicitly, indirect detection 
experiments. Three effects are a priori relevant here. First the large scale magnetic field of the 
Galaxy may prevent the DM particles penetrating the Galaxy plane if their gyromagnetic radius is 
smaller than the height of the disc. Concretely, for roughly 

a charged DM particle from the halo can not penetrate the Galactic disk [27]. Here B is the mean, 
large-scale magnetic field, and vdm the velocity of DM. A priori this may be counterbalanced by 
diffusion of DM by the small-scale, turbulent component of the Galactic field. The diffusion time 
scales like 

18 / Hd \ V GeV \ /300km/s\ V B \ 
Tdiff ^ W'^K ^ yrs . (40) 



with Hd the height of the Galactic disk. If we compare this with the propagation time scale 
is Hfi/vDM ~ 4 • 10^ yrs, we get that diffusion is effective compared to magnetic shielding only 
if K < 4 • 10~^^, a very small value which puts the model beyond the reach of direct detection 
experiment. One alternative is to assume that DM is present ah initio within the Galactic disk. 
However, this component of DM may be expelled from the Galactic disk by Fermi acceleration 
in supernovae shocks |33| |27] . The time scale for DM acceleration by supernovae shock waves is 
estimated to be of the order of Tacc ~ 10^ years |27) . Expulsion is efficient if either the time scale 
for energy loss, Tej,, or the diffusion time scale T^iff, is large compared to the acceleration time. The 
former gives k < 3.4 • 10-^(mDA//GeV)i/2^ while the latter requires k > 10 ^^{rriDM/ GeV), using 
the fiducial values of the parameters in Eq.(40). Hence, at the end of the day, we must assume that 

is required in order to match the expected local energy density, poM ~ 0.3 GeV/cm^. A glance at 



Fig. 10 reveals that this condition only concerns candidates with niDM ^ 2 GeV. 



Summary of cosmological constraints for a massless hidden photon: Taken literally the constraints 
discussed above point towards freeze-in as the only possible testable (z.e possibly no screening by 
magnetic fields) or viable (compatible with, say, the ellipticity bound) scenario, except for values 
of m£,M larger than a few hundreds GeV, for which the reannihilation and hidden sector freeze- 
out regime scenarios are also possible. This also excludes the possibility that the model would be 
responsible for the DAMA or CoGeNT data (regardless of the exclusion set by XenonlOO). These 
require moM in the few GeV range, and a At ~ 10~^, values that can lead to the observed relic 
density only for a' ~ 10"'^, which is much larger than what is allowed by the ellipticity bound, 
a' < 8 • 10~^. In other words, taking into account the ellipticity bound, such candidates would have 
a relic density about two order of magnitudes larger than what is observed. 
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Ellipticity bound and depletion for a massive hidden photon: Note that the bounds of Ref. |10| \TT\ 
[T2j apply exclusively for a massless mediator. One simple way to relax the ellipticity bound and/or 
the impact of magnetic fields, and so expand the scope of our scenario, is to assume that the C/(l)' 
group is spontaneously broken in such a way that the hidden photon becomes a light 7',^^ To 
assume that the 7' is massive does not affect either the relic density analysis above as long as 
m.y' < rriDM, or the direct detection analysis above, as long as my < 1 MeV. But it does relax 
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largely the cosmological constraints just discussed. To show that, we adapt the analysis of |Tl] to 
the case of a massive hidden photon Our results are shown in Fig. 11 
infer that, for moM ~ 10 GeV, taking ~ 1 MeV allows to take a' 



As an example we may 
10"^, which puts such a 

candidate in the reannihilation regime for the production of the DM abundance (see Fig. |3]). Note 
that the red curve in Fig. 1 1 corresponds to a massless hidden photon and has to be compared with 
Fig. 1 of [TT]. The DM candidates below the red curve are excluded. Our bound is compatible, 
but slightly milder than in [llj . This, we believe, is because we have taken an average DM energy 
density pDM ~ 1 GeV/cm^ and velocity dispersion vq = 300 km/s, rather than profiles for these 
parameters as in [llj. However this approximation serves our purpose which is to illustrate how 
the constraint from ellipticity changes for a massive hidden photon. Notice also that the constraint 
from interaction with the magnetic fields in the galaxy and supernovae shocks quickly drops if the 



7' is massive. This stems from the fact that the correlation length m 



< 



1 MeV"^ of the 7' is 



much smaller than the typical particle separation length, 0(1 cm). 

To show how the mass of the 7' affects the direct detection analysis above, let us consider for 



the massive case the differential cross section for scattering on a nuclei, substituting in Eq. (34) 
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da 



dER 



oc 



1 



1 



da 



dER {ER + ml/lmAY' 



oc 



(42) 



For A =Xe and Er in the few keV range, the mass of the hidden photon becomes significant 
only for my ~ yJmxeER > 50 MeV. From Fig. [oj we see that a candidate with m^M ~ 10 GeV 



^^The fact that the U (1)' gauge group is spontaneously broken does not affect the stability of the lightest fermion 
of the hidden sector. But still, even if broken, it forbids a possible Yukawa coupling to a SM lepton doublet and a 
Higgs doublet, which otherwise would make the DM particle unstable. 

^^From the cross section for the scattering between DM particles interacting through a massive photon, 



da 



4a 



m-DM 



the rate for energy loss in the Born approximation is given by 



E = ^^fp'^'' I dvve' 



log [AmjjMV^ sm^{0/2) + 



4mf 



where the minimum scattering angle 9min is set by Debye screening in the hidden plasma The integral over 

DM velocity v is peaked at = wo, so that for a massless 7', this reduces to the result of |llj . with the familiar 
divergence from Rutherford scattering, E oc — log(sin^(&min/2). The divergence is cut-off for a massive 7', E oc 
— log(m^//i;m_DM) but the dependence is essentially logarithmic for small m^/, so the impact is moderate unless 
rriji > few keV (see Fig. 111. The turnover behaviour of the ellipticity bound is related to the change of sign of the 
log at for m-yi ~ moMVo- 

^*One caveat is that our calculations are done in the Born approximation. Born approximation is only reliable for 
"weak coupling", ie a'm^i /mnMv'^ < 1 (see |39p. a condition which does not hold over the whole parameter range we 
consider. In particular, for m^i — 1 MeV, our calculations are reliable for tudm ^ IGeV, while for m^/ = 10 MeV, 
weak coupling approximation is only valid for a more restricted range, 10 GeV < moM ^ 300 GeV. Although the 
general trend of giving a mass to the hidden photon should be clear from the figure, one should keep in mind that, 
in the strong coupling regime, the scattering amplitude may possibly get a resonant enhancement, leading for some 
parameters to more stringent constraints than those shown in the figure. 

^^See also [371 136| for recent work on long range forces in direct detection. 
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Figure 11: Ellipticity bounds (see text) in the a' — tyidm plane, for interactions through a massless 
hidden photon (lower hue), and m^/ = 1 keV, 10 keV, 100 keV, 1 MeV, 10 MeV, 100 MeV (bottom- 
up). 

and K ~ 10~^ is below the XenonlOO exclusion limit, but within reach of XenonlT. Hence, light 
mediators open the possibility to also probe the reannihilation regime of the model, while being 
consistent with laboratory and astrophysical constraints. 

Other constraints for massive hidden photons: There are numerous astrophysical, cosmological and 
laboratory constraints on massive hidden photons and their coupling to ordinary matter. Here we 
may refer to Fig. 10 of Ref. |35] (see also |30) ) where constraints on the mixing parameter (our e, their 
x) and hidden photon mass are summarized, "^^ The range of interest for our discussion corresponds 
to m^' > 100 keV. Below this value, there are strong constraints on the mixing parameter from 
constraints on energy loss in stars. For 1 MeV < m^/ < 100 MeV, there is a constraint on e < 
10^^ from the duration of the SN1987A neutrino burst, which is consistent with the values of 
K ~ ey^ a' /a ~ 10~^ and a' ~ 10~^ required in the reannihilation regime. Above m^i ~ 100 
MeV, the constraints are much milder, but, in our opinion, we may not really speak of a light 7' 



anymore] I If the hidden photon is in the MeV range, one should worry about constraints from 
nucleosynthesis, possibly ruining the production of light elements [6]. For the sake of illustration, let 
us consider again the reannihilation scenario for ttidm ~ 10 GeV, and so e ~ 10~®. To begin with, 
we may be confident that the Universe is radiation dominated at the time of BBN, as the energy 
density in the DM and in the hidden photons are comparable around the MeV scale. Furthermore, 
T' <^T (for instance T' ~ T/10 for rriDM ~ 10 GeV and a' ~ 10^^), so that a massive photon in the 
MeV mass range becomes non-relativistic way before the epoch of BBN. Then, for m^/ in the MeV 
range, the only relevant decay is into e"^ — e~ pairs, with a decay rate V ~ ae^m^i (see for instance 
^), the lifetime of a hidden photon is C'(10~^s), and so takes place well before nucleosynthesis and 

""in the same reference [35] one may find in Fig. 9 various constraints on milli-charged particles summarized. There 



'Lighter hidden photons m^i < 1 eV are also less constrained but they have Httle impact on the elHpticity bound. 
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thus is harmless [6j. One last potentially interesting consequence is for indirect detection. While 
annihilation in visible sector degrees of freedom is suppressed by the mixing angle in the case of 
massless hidden photons, in the massive case, DM-antiDM pairs may annihilate in a pair of on-shell 
massive hidden photons, which subsequently decay into electron/positron pairs with energy mDM/2 
(the secluded scenario |6j). In the reannihilation and annihilation in the hidden sector regimes, the 
annihilation cross section is canonical av ~ 10~^^ cm^-s~^. For heavy DM candidates, a substantial 
boost is necessary to lead to observable consequences, but for lighter ones, moM ^ few GeV, there 
may be strong constraints, for instance from synchrotron radiation (see for instance |41| and |42)). 

Hidden atoms: Another easy way to evade the cosmological and astrophysical bounds for a massless 
hidden photon is to consider that there are more than one species of particle (for instance two, say 
e'l 2) in the hidden sector, which may combine into neutral states after freeze-out. Both particles 
would be stable because QED conserves flavour. Positronium-like bound states would still annihi- 
late, but the bound state e'^ 2 ~ ^2 1 ones would still be stable. If the binding energy is in the keV 
range, direct detection would probe the constituents ei,2 particles, but we have checked that the 
value of a' required to have recombination of 61^2 in the early universe is much larger than the ones 
required to have the right relic abundance, so we do not envisage this option as a viable scenario 
within our (minimal) framework. 



6 Creating a DM hidden sector through the Higgs portal 

Besides the kinetic mixing portal, there is another way to connect two sectors charged under different 
gauge groups: through the Higgs portal. Recently it has been considered in many different contexts. 
It is interesting to see how, in the same way as above, a DM hidden sector could be generated from 
the SM sector through this portal. In particular, we may wonder whether there is some universality 
in the structure we have seen emerging in the previous sections (the Mesa diagram)? This will turn 
out to be the case, even though in details there are some new phenomena related to the fact that 
the mediator, here the Higgs, may be heavier than the DM candidate. 

To this end we consider the simple case where the DM consists of a scalar charged under an extra 
[/(I)', in which case it can couple to the SM through the Higgs portal interaction, with coupling 

C 3 D'^^'^D'^<j) - Xm^^'^H^H - ^il(t)cp^ - A^(#"^)2 - ^i^H^H - \{H^ Hf , (43) 

with H = {H^ , {h+v) / We assume that the scalar has no vevso that the U{1)' gauge symmetry 
is unbroken and the scalar is stable. Neglecting the possibility of a kinetic mixing portalj^ the 

^■^In this scenario, we could consider the possibiUty to have simultaneously the kinetic mixing and the Higgs portals. 
However, in practice, unless the parameters of both portals are finely adjusted, only one portal should be relevant 
in a given regime. As our purpose here is more to confront a massless (or essentially massless) mediator to a heavy 
one, we do not consider further this possibility. 
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processes which can create the (p from the SM sector are 
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where i = 1, ...,n denotes the number of degrees of freedom of (j). In the following we will consider 



the Boltzmann equation for the charged scalar of Eq. (43), n = 2. Since the Boltzmann equation is 
the same for all real components of (j), the DM relic density for another value of n can be obtained 
by multiplying the DM relic density by n/2 (up to a small logarithmic dependence in the number 
of degrees of freedom). 

The main differences between the kinetic mixing and Higgs portal cases come from the fact that 



in the later case the mediator is quite massive, uih > 114 A GeV,^^ The consequences are two-folds, 
(a) DM can be also created by the decay of the mediator and (b) the various production channels 
are suppressed at small temperatures, either Boltzmann suppressed if the Higgs boson is real in the 
process, or by the mass of the Higgs boson at the fourth power if the Higgs boson is virtual. As we 
will see this implies important differences between the kinetic and Higgs portals but, still, in both 
cases there is a characteristic "Mesa" shape phase diagram. 

In the same way as for the kinetic mixing portal, for all phases but for freeze-in, it is required to 
determine p' / p as a function of T. For rriuM > Trih/'^ the DM is exclusively produced by scattering 



processes and, since it is produced in pairs, Eq. (10) applies. In this case p'/p scales in the same 



way as for the kinetic mixing, i.e. as ~ 1/T (~ const) for T > uidm {T < moM)- This can 



be seen in Fig. 12 which gives p'/p as a function of T for various values of itidm and A^- For 



TTT-DM < '"T'/i/2, on top of the energy transferred through the various scattering processes, given by 
Eq. (10), there is also a contribution from the Higgs boson decay, which turns out to dominate. The 
energy transferred by a decay is given by 

^ + AHp' = / f{p)E,r,{E,)d'p = ^T^mlTK^C^) , (49) 



or in terms of p'/p and T 

difilp) 1 



where Th{Ep) and Th stand for the decay width to DM pairs of a Higgs boson of energy Ep and at 
rest, respectively. This gives p'/p ~ 1/T^ for T > For T < mh («-e. typically for T < 7n,/j/6), the 
DM pair creation, and therefore the increase of p'/p, basically stop as the Higgs mediated channels 



^In the following for all numerical results we take the value = 120 GeV. 
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Figure 12: Evolution of the ratio of the visible and hidden sectors energy densities, for a range of 
connector parameter, Am = 10^^'^^'^'^'^^'^^'' (from up to down), and for various DM masses. 



become very suppressed. This is distinct from the kinetic mixing case where, because the mediator 
is massless, there is no such cut-off. These properties can be seen in Fig. 12P'^ The thermalization 



condition for the connector, Tconnect > H, which gives T'/T ~ 1, taken at T = ttidm, translates 
into Am > 5 • 10"^ for tjidj^j below ~ itiii/2 and Am > 6.3 • 10^'^ for uidm = ITeV. If Am is small 
enough for the coimector not to thermalize, p' / p reaches a plateau at T ~ Max[m£)M,i^h/Q] with 
value, p'/p ~ 7 • W'^X'^mpi/mDM {p'/p ^ W^Xl^mpi/mDM) for uidm < •mh/2 {moM > ruh)- 
In the intermediate ruh regime a rough estimate is p'/p ~ 10^^ X^mpi / mw . 

From the knowledge of (p' /p){T), one can integrate the DM number density Boltzmann equation, 
which takes the same form as the kinetic mixing one, Eq. (|6j), but since the decay term dominates for 
T^DM < "i/i/2 it is useful to write it down explicitly, subtracting it from the scattering contribution 



^"""^ = - (l - ^) + E (l - ^) +lHs{l- ^) , (51) 
where Y = n^j s = n(f,*/s stands for the number density of (p or (j)* particles (Yd A/ = 25^), with 
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Note that at T >> rrih the energy transfer is nevertheless dominated by the scattering terms because, as shown 



above, they give p'/p ~ 1/T. This explains the change of behaviour around T ~ 1 TeV observed in Fig. 12 
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Figure 13: Higgs portal phase diagrams for moM = ^e, 0.1 GeV, 10 GeV, 1 TeV (same label 
definitions than in Fig. |3j and the dashed blue line stands for the transition between IIIA and IIIB, 
as explained in the text). 



^connect 



Tig r/j(r) the decay reaction density of the Higgs boson to (pcj)*, and 



TUT) = r{h 



Ki{mh/T) 
^K2{mh/T\ 



The phase diagram, Fig. 13 turns out to have the same characteristic "Mesa" shape as for the kinetic 
mixing portal, despite the fact that the mediator is massive. It also displays 4 phases: freeze-in, 
reannihilation, freeze-out in the hidden sector and connector interaction freeze-out. Similarly, as 



for the kinetic mixing case, one gets a truncated volcano shape in Fig. 14 where Y is displayed as 



a function of Am for various values of a' and rriDM- There are nevertheless important differences 
one observes by looking closer at the way these diagrams are obtained. 
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Figure 14: DM relic abundance Ydm as a function of the connector parameter Am for different 
DM masses rriDM and values of the hidden sector interaction, log^o(a'/a) = 1, —1, —3, —5, —7, —9, 
bottom- up (the last two hues are the same for rriDM = 10 GeV, as well as for rriDM = 1 TeV). 



6.1 Phase I: the freeze-in regime 

The freeze-in production of a pair of scalar DM particles through the Higgs portal has been already 
considered in details in Ref. |l3] for = moM « "T-h, see also Ref. |44| . In this case the decay 
process dominates the DM production because it involves a smaller number of couplings and is not 
more Boltzmann suppressed than these scattering processes (unlike in the freeze-out case). The 
scatterings are responsible for less than 1% of the total DM production. Since the decay production 
is infrared dominated (until the number of Higgs particles becomes Boltzmann suppressed) one finds 
that the number of DM particles produced is simply the number of Higgs bosons times the decay 
rate times the Hubble timePj times a constant c of order unity 



nl^T{h DMDM) 



Y _ 

sH 



T=mu 



(52) 



c is independent of the decay, c = 37r/(2ii'2(l)) — 2.9. Here too c is bigger than unity because the 
maximum value of the production rate n'f^T{h — t- DMDM) occurs at T ~ m/i/3.5 rather than T ~ 
m/i, see Fig. 2 of Ref. |43| . Note that the DM abundance Y generated by the decay is independent 
of m£)M (modulo channel threshold effects) since its production stops at a temperature within the 



^^Notice that, in the freeze-in regime, the production of DM is indeed proportional to the decay rate of the Higgs 
boson into DM particles, and not to the branching ratio, because the Higgs boson number density is the equilibrium 
one, no matter the Higgs decays to SM particles are dominant. What is relevant is the competition between DM 
production from a equilibrium number density of Higgs boson and the expansion rate. 
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Figure 15: Higgs portal parameter required to get the observed DM relic density through freeze-in 
(a' = 0). 

range mh ^ T > rriDM, and because the Higgs decay width is independent of rriuM (to lowest order 



V /mh, 



the 
has 



in mjjj^/mj^). Therefore, taking into account the dependence r(/i — )• DM DM) oc A, 
DM relic density scales as m^A^/m^ and, the larger the Higgs boson mass is, the larger m^A, 
to be to reproduce the observed relic density. This parametric dependence can be seen in Fig. 15 
which gives the value of A^ necessary to get the observed relic density as a function of rriDM (see 
also Fig. 14). Note that in this case the DM particles are produced when they are relativistic, since 
m/i/2 > rriDM, but as well-known, as long as moM is above the ~ keV scale, this is still compatible 
with the constraints on the fact that DM should be cold. 

In the opposite regime, nioM > nT'h/'^, the Higgs boson decay production ceases to contribute 
and the scattering processes become the relevant ones, producing DM down to T ~ rriDM- This 
gives in a similar way 



Y 



^connect 



sH 



(53) 



with c of order unity. For the WW, ZZ, hh channels we get c = 2.5, 2.5, 2.4 respectively. This leads 
to a relic density, i^oM oc A^, independent of ttidm- As a result, as can be seen in Fig. 15 the 
value of Am which gives the observed relic density. 



A. 



7.6-10 



-12 



(54) 



is also independent of moM in this regime. The WW process is dominant but is not the only relevant 
one. For example for = 1 TeV the WW, ZZ,hh and tt processes contribute respectively to 

49%, 25%, 25% and 0.5% of the relic density. 



6.2 Phase II, III and IV: the reannihilation and freeze-out regimes 

In the case of the kinetic mixing we have seen that, if the hidden sector thermalizes, but not the 
connector, the DM abundance freezes according to the reannihilation regime, where both the source 
term and the hidden sector thermalization process are active down to the freezing temperature. In 
this explained above, Y first follows the thermal density Yeq{T'), then follows the quasi 

static number density Yqse until it intercepts Ycrit where it freezes. Y follows Yqse because Yeq{T') 
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Region.III-A(ffJ[3M= lOGeV, A,„=10"'', = 10"'') Region.II (m [5^= lOGeV, A„,=10"", tr' = 10"") 




0.01 0.1 1 10 100 0.01 0.1 1 10 100 



z='«dm/T z=W1dm/T 

Figure 16: Examples of evolution of the DM number density K as a function of z = uidm/T. The 
first example corresponds to a case of hidden sector interaction freeze-out without thermalization 
of the connector (region IIIA), obtained with uidm = 10 GeV, Am = 10~^ and a' = 10~^. The 
second example corresponds to a reannihilation case (region II), obtained with = 10 GeV, 

Am = 10~^^ and a' = 10"^''^. Also shown are Yqse (blue), Ycrit (purple) and Yf.q{T') (green). 



intercepts necessarily Yqse before it intercepts Ycrit {Yf>q{T') gets suppressed at T' < moM, that is 
to say before Yqse, since the later {i.e. the connector) is suppressed at T < itidm)- 

For the Higgs portal, if rriuM > the situation is the same. For example for raoM = 1 TeV 

in Fig. 13 we get a phase diagram which has the very same structure as for the kinetic mixing 
interaction. However if rriDM < nT-h/'^ the situation is different. In this case, instead of being 
suppressed at T = m£)M, as for the massless mediator case, the connector interaction (and hence 
Yqse) gets cut-off at T < mh/Q- As a result, unless T'/T is small, Yqse gets always suppressed 
before Yeq{T'), since the later gets suppressed only at T' < niuM- In other words the system passes 
through a period with nioM < T' < T < mh/Q, where the connector has already decoupled but not 
the hidden sector interaction. An example of evolution of Y where this holds is shown in the left 
panel of Fig. 16 There one sees that, as the hidden sector thermalizes, Yeq{T') gets larger than Yqse 
(as in the massless case) and then stays so, i.e. Y follows Y^qiT') until it intercepts Ycrit without 
ever following Yqse- In this case there is no period of reannihilation (where Y follows Yqse), even 
though T'/T < 1, but a hidden sector interaction freezes-out. This freeze-out is standard in the 
sense that it occurs when the connector has already decoupled {i.e. when T'/T is already constant). 
But still, the final relic abundance depends on the size of the connector because, prior to freeze- 
out, the larger the connector is, the larger T'/T is, and so the larger is i^DM- This regime has no 
equivalent in the massless case. It is a regime which comes between the reannihilation regime and 
the hidden sector freeze-out regime with thermalization of the connector. In the later case, the relic 
density is independent of the size of the connector (since T'/T ~ 1). 

Practically it is clear that without thermalization of the connector, the lighter is uidm, the 
easier the system can evolve through a situation where muM < T' < T < rrih/G, hence the larger 
is the region of parameter space where this freeze-out phase occurs, and conversely the smaller is 
the reannihilation region. For instance for moM below ~ 1 GeV there is basically no more relevant 
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reannihilation phase region, i.e. the reannihilation occurs with values of Am too tiny to give a large 
enough relic density. Similarly for rriuM ^ "^/i/6, there is, as with the kinetic mixing interaction, 
no sizable region of parameter space where the hidden sector interaction would freeze-out without 
thermalization of the connector. As for the intermediate DM mass values, they allow both regions. 
An example of reannihilation relic density evolution for an intermediate mass value, moM = 10 GeV, 



is given on the right panel of Fig, 16 Let us now discuss each regime with equations, for the case 



Reannihilation : region II in Fig. 13 for the rriDM = 10 GeV panel. Starting from a freeze-in regime. 



and increasing Am, the hidden sector gets sufficiently populated to thermalize For values of 
parameters where thermalization is hardly obtained, eventually T'/T will still be sufficiently small 
for Yeq{T') to be suppressed before Yqse- As a result, after having followed YggiT'), Y follows 



the Yqse line, until freeze-out, see right panel of Fig. (16). In this case reannihilation proceeds as 
in section 3.2 above, except that one has to replace everywhere the {(Jconnectv) scattering term by 
the decay term {'^connect = {crconnectv)nlg(T) Th(T)n'^g). In particular Ycrit(T) = H/{aHsv)s is 
unchanged but Yqse for a decay is given by 



YlsEiT) = Y,„tTH{T)^Y,,{T) = ^I^n%{T) . (55) 



At freeze-out, where Y = YQSE{Tf) = Ycrit{Tf) one gets therefore, in the same way as Eqs. (26)- 



(27), 



^n^q{Tf)ThiTf){aHsv) = H{Tf) . (56) 

This equation can also be found in Ref. H]. Taking into account the fact that Ty < rrih (when the 
connector source term gets suppressed), i.e. n'^g(T) = (^^^^)^/^exp[— m/i/T] and Th(T) ~ T^, from 



Eq. (56), we get 



Xf ~ log 0.0229c c + 2 ^TF^^b^ 

+ - log[log[0.0229c(c + 2) ^— ]] , (57) 

^ 9* m^V 

From this value of x/ one gets therefore Y{Tf) by substituting it in 

YiTf) = Y^ruiTf) = —V = — 7 ^ (58) 

{(THSv)s T=Tf 2, 2ti^ g^s {crHsvmDMmpi 

c is a numerical constant of order unity. One gets the numerically exact result for c between 0.4 
(for large values of a') and 1 (for small values of a'), related to the fact that Y freezes-out in reality 
a little bit later after that Yqse = Ycrit)- In this regime we find that Xf varies between ~ 1 and 
~ 5. As with the kinetic mixing interaction in the reannihilation regime, the final relic density is 
inversely proportional to (ansv) but Xf depends on both interactions, so that it is smaller than the 
usual ~ 20 value. 

The transition line between the freeze-in and reannihilation regimes is given by the condition 



H 



~ 1 (59) 

T=mr)M 



^Unless a' is so small that it never thermalizes even if T'/T ~ 1. 
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which using Eq. (52) can be rewritten as 



H 



T=mDM 



r,(rK,(r)|^^^^^) c^i (60) 



27 



with T = Max[m£)Af,~ m/i/6] 

Note finally that here too, the DM particles are produced relativistically from the decay of the 
Higgs boson, but, unlike in the freeze-in regime, quickly they become non-relativistic from the fact 
that they are in kinetic equilibrium with the 7' down to a T' temperature below muM- 



Hidden sector freeze-out without thermalization of the connector : region IIIA in Fig. 13 Starting 



from a reannihilation or freeze-in regime situation, see Fig. 13 if one further increases \rn^ some 
point T' gets sufficiently high for the hidden sector scattering not to be much Boltzmann suppressed 
when the connector T ~ mh/6 decouples, even though T'/T < 1. In this case the hidden sector 
freezes-out when the connector has already decoupled for a while. Therefore to get the freeze-out 
temperature as well as the relic density, one just needs to calculate the final value of p' / p which is 



reached at T ~ rrih, see Fig. 12 and is approximately given by 

p' Th{T)-mhn%Q{T) 



Hp 



Tr-.^ (61) 



which gives ^ = T'/T = {gf^'' / g^^Y^^ip' / pY/"^ , and to plug this value of ^ in Eqs. (|32(-(|33[. This 
gives a relic density which scales as m^j^^Am^/a'^- The separation line between the reannihilation 
regime and this regime is given by the condition 

Yqse = Yeg{T') (62) 

at freeze-out where Yf>q(T') = Ycrit- Numerically it corresponds to Am ~ 3 • W~^{m£)M/'mh)'^ 



(neglecting a logarithmic dependence on a'). This condition, together with the condition of Eq. (60), 
leaves relatively little parameter space for the reannihilation regime. For example for rriDM = 
10 GeV, reannihilation occurs within 2-10"^^ < Am ^ 2-10"^'^ (imposing the relic density condition). 
Practically this means that apart for the region 1 GeV < < mij?,^ the reannihilation regime is 

not relevant to give a large enough relic density and this regime fills all the parameter space where 
the hidden sector interaction thermalizes but not the connector. Note also that the value of we 
get in this case lies between ~ 20 (for T'/T close to 1) to ~ 5 (for lower values of T'/T). 



Hidden sector freeze-out with thermalization of the connector : region HIB in Fig. 13 If one in- 
creases further Am, at some point, the connector interaction tliermalizes, wlien ^connect — 

H at 

T > Max[mj:)M-, ~ mh\^ and one enters in the standard hidden sector interaction freeze-out regime. 
In this case ^dm scales as ~ ^/{cTannih'v) oc m^j^^/a'^ and ~ 20. It is independent of Am 
since T'/T ~ 1, at least in first approximation. In second approximation one can mention never- 
theless that the relic abundance might be not totally independent of the size of the connector for 
T^DM « nT'h- This stems from the fact that, as the communication between both sectors stops at 
T ~ mh/6, the hidden sector temperature cannot track the visible sector one when T gets reheated 



^^One can also mention that for intermediate DM mass values, at the transition between freeze- in and reannihilation, 
the relic density undergoes a "kick", see Fig. |13| a feature which is not found in the kinetic mixing case. It can be 
shown that it is due to the fact that it corresponds to a transition from a situation where DM is relativistic (in the 
freeze-in regime, with a kinetic energy of order ~ mh/2) to a situation where DM, as soon as it has thermalize, looses 
a lot of kinetic energy (down to a kinetic energy of order T'), leading to a enhanced hidden sector annihilation cross 
section, since {(Jhsv) oc 1/s, hence a reduced DM abundance. 
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due to the decoupling of a SM species (but this is a rather small effect, as can be seen in Fig. 13 
for rriDM = and Am ~ 10~^). 



Connector freeze-out : region IV in Fig. 13 Finally if the connector not only thermalizes but 
also leads to a reaction rate larger than the one of the hidden sector, then one enters the standard 
connector freeze-out regime. For a relic density of order of the observed one it goes through the usual 
freeze-out of the annihilation process with H-dm ~ / {o'connectv) ■ In particular, for mohi < mh/'^i as 
is of order 20 the Higgs decay-inverse decay process can be totally neglected at freeze-out since it is 
obviously more Boltzmann suppressed than the scattering process. This regime has been abundantly 
studied in the literature for the Higgs portal, for the case of a real scalar (n=l) |45| Il6l 14 7| . as 
well as for n = 2 |48j . For uidm ^ ^hf^ it gives a relic density which scales as m^/Am (for 
rrih ~ 120 GeV) whereas for tudm « f^h it scales as mf^/{m'jX'^) with / the closest in mass 
fermion (with ruf < uidm) (and for m£)]\j » rrih it scales as mj-)j^/X'^). 

Note nevertheless that, for larger values of i^DM, freeze-out in this regime occurs for smaller 
values oi Xf, in which case the Higgs boson decay/inverse decay processes dominate the freeze- 
out process. This explains why along this regime the transition line between region HI and IV, 



corresponding to the condition Tamnect = ^HS: bends towards smaller values of Am, see Fig. 13 for 
the rriDM = 10 GeV example. 

Note also that for a' ~ 0, due to this bending behaviour, starting from Am — and increasing 
Am, one lays first in the freeze-in regime, then in the connector freeze-out regime dominated by 
the decay {xf ~ 1), then in the connector freeze-out regime dominated by the scattering processes. 
As a result in Fig. [14] the volcano is not perfectly conic as for the kinetic mixing case but remains 
truncated even for negligible values of a'. This behaviour is nevertheless irrelevant for cases which 
give the observed relic density. 

In the same vein, note also that, as for the kinetic mixing case, all I-II-III-IV phase transition 
lines meet at a single point. The meeting point corresponds to the situation where the hidden sector 
interaction, as well as the connector interaction, decouples as soon as they thermalize. Within the 
mass range we consider, me < niDM < W TeV, this situation always gives a too large relic density. 

Finally let us mention that for rrioM ^ 'nih/2 the III to IV transition line can be significantly 
modified because in this case the connector scattering process, which determines the DM freeze-out, 
is resonantly enhanced. 



6.3 Testing the Higgs portal phase diagram? 

To test the freeze-in scenario, as well as the reannihilation scenario, from direct detection is much 
more difficult for the Higgs portal than for kinetic mixing, because the former does not lead to 
a collinear enhancement of the direct detection cross section. For instance, the sensitivity 

reach from non-observation of a 100 GeV scalar at XenonlOO (XenonlT) is Am ^ 10^^ (resp. 
2 • 10~^), which is several orders of magnitudes larger than the characteristic couplings required by 
the freeze-in or reannihilation scenarios. This upper bound nevertheless constraint the Higgs portal 
interaction freeze-out regime which requires typically Am ~ 10~^ (Am ~ 10~^) for rriDM ~ 100 GeV 
(resp. 10 GeV) (for rrifi = 120 GeV, see e.g. Ref. |461 147) ) . as well as part of the hidden sector 
freeze-out regime. Similarly the freeze-in and reannihilation regimes cannot be tested either from 
collider experiments but possibilities of tests at colliders for large value of Am do exist, in particular 
from the invisible decay width of the Higgs boson, see Ref. j^9t 147) . In other words even for DM 
masses above the GeV scale the phase diagram is much less constrained experimentally than in the 
kinetic mixing case. Only the connector freeze-out regime and perhaps a part of the hidden sector 
interaction freeze-out regime could be tested within a not too far future. 
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7 Summary 



We have considered the simple possibihty that, along the thermal history of the Universe, DM 
could have been created out of SM particles through a portal connecting the visible sector (SM) to 
a dark sector. This framework is particularly predictive as it is induced by the SM particles, that is 
to say particles whose properties (mass and interactions) and thermal number density are known. 
As a result, the energy transferred from the visible sector to the hidden sector and the DM relic 
density produced, depend, besides (obviously) the DM mass, only on the strength of the connector 
interaction and of the interaction(s) which may thermalize the hidden sector. 

Whether we consider the gauge kinetic mixing portal, which involves a massless (or massive) 
mediator, or the Higgs portal, whose mediator is the Higgs doublet, we have found that there 
are 4 distinct ways of producing DM, leading to a phase diagram with the characteristic shape of a 
"Mesa". In two of the phases, the DM relic density depends exclusively on the connector interaction, 
whereas in the two others it depends essentially on the strength of the hidden sector interaction 
(provided the connector interaction is within a certain range). 

It is useful to summarize the parametric dependence of the relic density for the various phases, 
which holds for a connector lighter than the DM particle (both in the case of kinetic mixing and 
Higgs portals): 



I : Freeze-in: Y 



{n^eq^f{cFconnectv) 



sH 



with Xf ~ — ''^DM /ggs 

T=Max[msM,mDM] Max[msM , W-DmJ 



II : Reannihilation: Y — — ^ with Xf ~ log(cJe/-fv)?7i^ (64) 

{aHsv)fh^ 

III : Hidden sector freezeout: Y ~ \ ~ 2 "^^^^ -^Z ~ ^og{aHsv)rh'^ (65) 



{(THSv)rn? 
Xf 

connect'^)'^ 



IV: Connector freezeout: Y^- — with x f ^ log{aconnectv)rh'^ (66) 



with m? = mnMfnpi- If none of the interactions thermalize, the freeze-in regime holds. If the hidden 
sector thermalizes but the connector does not, reannihilation occurs. In this regime both interactions 
freeze at the same time after going through a period of quasi static equilibrium evolution. We have 
shown that in this case the relic abundance is inversely proportional to the {(Thsv) hidden sector 
cross section, but that the freeze-out temperature Xf is determined by an effective cross section 
which is the geometric mean of both the hidden sector and connector cross sections (aeffv) = 
\/ {'^HSv){'yconnectv)- If instead, both interactions thermalize one gets an ordinary hidden sector or 
connector interaction freeze-out regime, depending on which of the two interactions is the fastest. 
Finally if the hidden sector interaction is so small that it never thermalizes, even if the connector 
thermalizes, one also lies in the connector freeze-out phase. We have also shown that there are 
no other possible regimes and that, for a non negligible hidden sector interaction, the four regimes 
necessarily follow each other in this order. 

In the opposite case where the mediator has a mass larger than the DM, as applies for the Higgs 



portal for rriDM ^ 'm-h/'^ '^^ one gets the same pattern except that there is an additional decay 
contribution and that the connector interaction is cut-off at the mediator mass. As a result the 
decay dominates both the freeze-in and reannihilation regimes, and one has also the possibility that 
the hidden sector interaction freezes out without thermalization of the connector. Correspondingly 

^*Or for the kinetic mixing portal with moM < m^i /2, a possibility we didn't consider here but which would give 
a same pattern. 
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one has: 



I : Preeze-in: Y ~ 



sH 



with Xf ~ — (67) 

T=mh rUh 



2 3/2 

II : Reannihilation: Y ^ ^- — ^ with Xf ~ log[(aHSv)Th ^L^^ ] (68) 

IIIA : HS freeze-out (no connector thermal.): Y ~ _ „ with Xf ^ flog(aHSv)fh^(69) 

{aHSv)m^ 



IIIB : HS freeze-out (connector thermal.): Y ^ -. — — ^ with xt ^ \og{aHSv)rh'^ (70) 

{aHsv)m^ 



IV: Connector freeze-out: Y — ^ with Xf ^\og{accm,nectv)'m^ (71) 

with the Higgs boson decay width to DM particles, and ^ = T'/T = {g^^ /g^^)^^^{p' /py^"^- The 
reannihilation, as well as the hidden sector freeze-out regime (without connector thermalization) 
give a relic density that depends mainly on the size of one interaction, but also logarithmically on 
the other one, unlike the other regimes for which the relic density depends only on the strength of 
one interaction. For a heavy mediator, the reannihilation regime is relevant only for a relatively 
narrow intermediate DM mass region below the mediator mass, 1 GeV < uidm ^ TUh/S. Therefore, 
apart for this region, if the hidden sector interaction thermalizes but not the connector, hidden 
sector freeze-out always occurs. This is different from the low mediator mass regime where in this 
case {i.e. if T' < T) reannihilation always occurs. 

One should also emphasize the fact that one also gets a characteristic Mesa phase diagram in 
the more complex situation in which there are both a massless and a massive mediator. This applies 
in particular in the kinetic portal scenario, in the intermediate mass range 1 GeV< uidm ^ 
see Appendix D. 

To discuss the kinetic mixing portal case, we have considered a particularly simple model where 
the hidden sector consists of a single particle charged under an unbroken U^l)' gauge group, depend- 
ing only on m£,M- on the kinetic mixing coupling and on the extra gauge coupling. In addition to 
the relic density phase diagram, we have discussed in detail the phenomenology it may imply. Direct 
detection is particularly interesting for this model because it involves a massless or hght mediator, 
which implies that the cross section is inversely proportional to the recoil energy squared and is 
therefore strongly enhanced at low recoil energies. We have shown that the latest XenonlOO data 
exclude the regimes of freeze-out (both in the hidden (III) and visible regimes (IV)) for candidates 
with niDM >~ 5 GeV (see Fig.9). The same data also probe (i.e. exclude) a large fraction of the 
parameter space corresponding to the reannihilation regime (II), but does not constrain freeze-in 
(regime I) so that, for the time being, the latter is allowed for any DM mass. However, although 
the kinetic mixing parameter required to produce the observed DM abundance through freeze-in 
is very small, even this regime could be tested by the future Xenon-IT experiment. For one year 
(four years) of exposure, we have found that candidates in the range 50 GeV < ttidm ^ 140 GeV 
(50 GeV < niDM < 600 GeV) could be tested. As for the DAMA and CoGeNT data, they can 
be accounted for by considering values of the parameters which lay in the reannihilation regime, 
but which are excluded by the current Xenon-100 data, see Fig.8. The characteristic l/E^ recoil 
energy spectrum allows to distinguish this model from more standard DM models which, involving 
a mediator beyond the ~ MeV scale, predict a cross section on nucleon that is independent of Er- 

Other important constraints for the kinetic mixing setup are cosmological, related to the long 
range interaction caused by the hght mediator. The most relevant constraint rests on the ellipticity 
of galaxies, which puts an upper bound on a' . There is also another stringent upper bound on a' , 
which is related to the galactic magnetic fields and which one must satisfy if one wants DM to be 
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present at the Sun location, so as to be probed using direct detection experiments. Both constraints 
turn out to be fully compatible with the freeze-in regime, but not with the other regimes, except for 
a DM mass larger than a few hundreds GeV. We have shown how these constraints get considerably 
relaxed if instead of considering a masslcss mediator we make it slightly massive (still lighter than 
the DM mass scale, so that the relic densities we have obtained are unaffected). 

The model we have considered for the Higgs portal consists of a scalar charged under a ^7(1)' 
gauge group, and depends only on moM, the U (1)' gauge coupling and the Higgs portal interaction. 
Here the 4 phases arc all allowed by experimental data, even for DM masses above the GcV scale. 
Direct detection together with collider experiments could potentially cover the connector freeze-out 
phase and, optimistically, a part of the hidden sector freeze-out regime. The phase diagram we have 
derived can be consequently considered as the analytic prolongation (towards smaller Higgs portal 
coupling values) of the Higgs portal freeze-out regime widely considered in the literature. 

A priori such an analytic prolongation could be obtained with the same characteristic mesa 
shape for any DM candidate which belongs to a SM gauge singlet hidden sector, with (a) a visible 
sector which consists of particles that are in thermal equilibrium when they produce the DM (as 
in the case we considered here, where in the visible sector there are only SM particles), and (b) 
a primordial energy density in the hidden sector which is secondary. In this sense our results are 
representative of a large class of new DM models. 
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Appendix A Relevant cross sections for the kinetic mixing portal 



DM pair creation from SM fermions: 

This is the dominant process for DM creation through the kinetic mixing portal. It takes place 
through a photon or a Z in the s-channel: 



- Amis - AmliS + IGrnlml, Ami + Ami, 
(1 + 1 ^-1 + f ^ (72) 



X 



cos 6/|y ' S — m\ ~ imz^z 

+ 4cosg^(.-m|)' + r.|r| [^^^ - ^^^"^/^'^ + ) 

- Amis- Ami, s + lQmlmi, 
+{9v + 9a){s^ + ^sml - 8ml,mj + ^ ^ ^)]}, 

where e = e cos 9w /Vl — , and gv, gA, ^'f are the V, A components and color factor of the SM 
fermions with / = e, fj,,T,u,d, s,c,b,t. For moM » mz the 7, Z and 7 — Z interference terms 
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involving charged / particles, and the / = z^e./^.r channels, contribute along the 8 
proportions respectively. 

Z-decay in DM pairs: 

A Z on-shell may decay into DM pairs, with a width given by: 



-3.0 : 0.6 



^ ) = ^— mzWl - ^(1 + ^ 



rrir 



rrir 



(73) 



DM/SM fermion scattering: 

In practice this t-channel process, which is dominated by photon exchange (Z contribution not 
shown), is negligible and is given for reference (see the discussion in Section 2): 



aUe' ^ fe') = ^ 



irqjqlaa'i^ 

-,2 /y-n'^ \ I 2 



^ s|Picm(m2,m2,)|2 y_4|Pi,^(m2,^2^)|2 

with |-Picm("ii, ?Ti2)P = (s — mf — 1712)'^/ (4s) — mlm^/s. 
DM pair annihilation into hidden photons: 



(s-m^-m2,)2 s i 



t2 



dt , (74) 



a[e e — 7 7 j 



47r(a') 



f\2 



2s^ + 8m^,s - 16m^, 
(s — Am?^,)s 



tanh 



-1 



4m2, 



s + 4m^, 
s{s — ^m?i) 



(75) 



Appendix B Reaction densities 

In Eqs. ([6| and (51) the reaction densities (and accordingly the {av)^s) are defined as 

We have defined dp = (fp/{{2Tr)^2E). Here f^'^ = (e'^'/^zbl)-^ ~ e"-^'/"^ is the Maxwell-Boltzmann 
energy distribution, |A^p is the amplitude squared summed over initial and final spins (with no 
averaging), Smin = max[(ma + rrib)'^, (mi + 1112)'^], and the reduced cross section is defined by 



a{ab^ 12) 



9agb 2[{s-ml-mlf -Amlml] , . 
Cab s 



(77) 



with a the particle physics cross section of Eqs. (72), g^fi the number of degrees of freedom of the 



particles a,h and Cab a combinatorial factor equal to 2 (1) if a and b are identical (resp. different). 



As for the Hubble constant and the entropy density they are given hy H = 1.67\/ gt^^T^/mpi, with 
gt" = <?f *^ + 9^', and s = (27rV45)ff,,r3, with 5*. = <?ff + 9^s' ■ 
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Appendix C Energy transfer Boltzmann equation 



It is convenient to rewrite Eq. ([9| as 
dt 

where 1521 



/d'^pi d?p2 
J^^J^^h{Pl)f2{P2)vMol£{Pl,P2) , (78) 



£iPi,P2)= 2E^2E2VM I (/ ^2^% ■ g^\iMi2^3i\\27r)^6'~^\pi + p2 - P3 - P4)AEtrj , 

with vmoI the so-called M0ller velocity |52) . 

VMol = F/{EiE2) = yj[s - (mi + ms)^] [s - (mi - m2Y]/{2EiE2) . (79) 

In these equations /j is the momentum distribution of the initial "i" particle, p' is the pressure of 
the dark sector, gi is the number of effective degrees of freedom, /S.Etr is the energy transferred 
in the process and A^i2-h>3 4 the usual invariant matrix element. Note that as long as T' /T is not 
close to one, it is a good approximation to neglect the inverse processes because they transfer an 
energy to the hidden sector that is suppressed by a factor ^'^ (one factor of from the hidden 
sector particle number density and one of ^ from the energy transferred per scattering). We will 
nevertheless incorporate these processes to get results that are reliable for large ^ values too. 

Assuming a radiation dominated Universe in both sectors down to T ~ 1 eV, i.e. p = p/3, and 
adopting Maxwell-Boltzmann momentum distribution one therefore obtains 

^ + 4//,' = «» / llll^-"""--*'" (») 

For the e~^e~ — )• e'~^, e'~ pair production process the Boltzmann equation integration is relatively 
straightforward because the energy transferred, AEtr = Ei + E2, is independent of the momenta 
of final particles, so that £{pi,p2) is precisely cr(pi,p2) ■ AEtr- Since the cross section is Lorentz 
invariant it can be calculated simply in the center of mass frame, which, defining E± = Ei ziz E2 
and s = {Ei + £'2)^ - {pi +^2)^, gives 

dp' , .TT I f d^Pl d^P2 -Ei/T -E2/T -A ATT 

-^+^Hp = g^g^J——.e 1/ e a{pi,p2)vMoAEtr 

j dE+dE^ds ■ e-^+/'^a{s)FE+ 

j ds j dE+E+e-^+/^a{s)F j dE^ . (81) 



9192 
327r4 
5152 



327r4 

Integrating on the parameter space given by | I50] 

s > (mi + m2)^ 
E+>^s 



/ Ei>mi\ 




E2>m2 


H 


\i P1P2 1 - 





IeI-s<E^- E^'^^i^ < ^-^JE 



and since mi = m2 = m, one obtains the simplified form 

dp' . .rr f 5152 /■ , f t:^ -E^/T^^ 



dt 



+ ^Hp' = I ds I dE^E^e-^^/^^-^,/^sa{s)F 

5152 f , / 

/ ds—sTK2{ — )a{s)F 



327r' 



5i_52 
32 



g I rf, . ^(,)(, _ ^^2^sTK20^) , (82) 
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which agrees with the result of other works, such as Ref. |5lj . In terms of p'/p and the variable T 
this gives Eq. (10). 



Appendix D Interplay of decay and scattering processes for the 
kinetic mixing portal 

The interplay between decay and scatterings processes in the kinetic mixing portal is similar to that 
in the Higgs portal except for the important difference that, for the former, not all scatterings are 
mediated by the decaying particle (the Z in this case). Indeed, for kinetic mixing in addition to 
the Z e'e' decay and Z mediated // — )• e'e' processes, there is also a contribution from 7, which 
is not suppressed at T << mz- As a result, while for the Higgs portal the Higgs decay always 
dominates the DM production for m£,M < '^/i/2, for the kinetic mixing this is true only within 
the range ~ 1 GeV < moM < mz/2. This is due to the fact that, as we have seen before, DM 
production from the 7 contribution to // — )■ e'e' is infrared dominated, i.e. it is enhanced at low 
temperatures (and maximum at T ~ Max[mDM,Tnf]), while the production from the decay occurs 
at a temperature that is a few times smaller than mz (before it becomes Boltzmann suppressed). 
It is useful to give a few more details about the interplay between decay and scattering processes 
that takes place in this intermediate mass range, and this for each regime. 

Freeze-in : In this regime the decay channel {i.e. the resonant part of the scattering) dominates 
both the energy transfer and the DM production. The discussion is exactly the same as for the 



Higgs portal and we will not repeat it. In particular, just in the same way as in Eq. (52), one has 

nf r(Z ^ DMDM) 



Y 



sH 



(83) 

T=mz 



with c = 37r/(2/s:2(l)) ^ 2.9. 

Reannihilation : In this regime the situation is slightly more subtle than for the Higgs portal. In 
the case of the Higgs portal the QSE line is dominated by the decay contribution for raoM below 
mh/'2. In particular, just like the contribution from decay, the scattering process is cut off at low 
energies by the Higgs boson mass scale. Similarly, for kinetic mixing the Z decay is cut off by the 
Z mass scale but the 7 mediated // — )• e'e' process is cut off only at T ~ Max[m£)j\,/, mj]. As a 
result Yqse niay display two bumps, one from the decay at T mz and one from the scattering 



at T ~ m£)M (or at T ~ m-e if m^M < nT-e)- This is shown for two examples in Fig. 17 both 
for moM = 10 GeV. Therefore, to determine how reannihilation may occur, the key issue is to 
know when Yqse crosses Ycrit (at which point the relic density freezes). If this happens on the first 



bump, then the decay rate determines the relic density, as in the first panel of Fig. 17, and the decay 



reannihilation equations of Section 6 apply. Conversely, if it occurs on the second bump, as in the 



second panel of Fig. 17 then the 7 mediated process determines the relic density and the scattering 
reannihilation equations of Section 3 apply. Actually, since Yqse scales as 1/a' while Ycrit scales 
as 1/a'^, large values of a' imply that the dominant contribution is from scattering, whereas for 



small values of a' the dominant contribution is from Z decay. In Fig. 18 we show, in addition to 
the possibilities described in the moM = 10 GeV panel of Fig. [3j a line which delimitates these 
two reannihilation regimes, with HA (HB) reannihilation through Z decay processes (respectively 7 
mediated processes). 

Freeze-out regimes: For kinetic mixing the situation in the intermediate mass range is exactly 
similar to that of the Higgs portal, there is IIIA, IIIB and IV regimes. All cases the production is 
dominated by decay and, once both sectors have thermalized, there is standard freeze-out of the 
annihilation processes. If the case of hidden sector freeze-out, in which the connector does not 
thermalize, the transfer of energy is dominated by the decay channel. 
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Region.n-A(mnM= lOGeV, ;t=10""'\ ff'=10~''^) Region.II-B (m lOGeV, ^=10"*', o-' =10"^ *) 




z=n?DM/T Z=mDM/T 

Figure 17: Examples of evolutions of the DM number density as a function of z, in the reannihilation 
phase with uidm = 10 GeV. The first (second) panel corresponds to a case where Ycrit intercepts 
Yqse on the Z decay bump (7 mediated scattering bump). Both cases give the observed relic 
density. The first panel is obtained with k = lO"^*^'^ and a' = 10~^'^ whereas the second is 
obtained with k = 10~^-^ and a' = 10^^-^. 



phase diagram Log|(,[yDMj (wdm^IOGcV) 




-12 -10 -8 -6 -4 -2 



Figure 18: Phase diagram for the kinetic mixing portal and uidm = 10 GeV separating explicitly 
the reannihilation regimes dominated by the decay (IIA) and by the 7 mediated scattering (IIB). 
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